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Abstract. Attractiveness is a fundamental tool to study interacting particle systems 
and the basic coupling construction is a usual route to prove this property, as for instance 
in simple exclusion. The derived Markovian coupled process {Ct,Ct)t>o satisfies: 
(A) if ^0 < Co (coordinate-wise), then for all i > 0, < Ct a.s. 

In this paper, we consider generalized misanthrope models which are conservative par- 
ticle systems on Z"* such that, in each transition, k particles may jump from a site x to 
another site y, with fc > 1. These models include simple exclusion for which k = 1, but, 
beyond that value, the basic coupling construction is not possible and a more refined one is 
required. We give necessary and sufficient conditions on the rates to insure attractiveness; 
we construct a Markovian coupled process which both satisfies (A) and makes discrepan- 
cies between its two marginals non-increasing. We determine the extremal invariant and 
translation invariant probability measures under general irreducibility conditions. We ap- 
ply our results to examples including a two-species asymmetric exclusion process with 
charge conservation (for which fc < 2) which arises from a Solid-on-Solid interface dynam- 
ics, and a stick process (for which fc is unbounded) in correspondence with a generalized 
discrete Hammersley-Aldous-Diaconis model. We derive the hydrodynamic limit of these 
two one-dimensional models. 



1. Introduction 

Attractiveness is a fundamental property to study interacting particle systems: combined 
with coupling, it is a key tool to determine the set of invariant (and translation invariant) 
probability measures for the dynamics; it is also essential in obtaining hydrodynamic limits 
for conservative systems under Euler scaling (see [18], [IS]). Most widely known examples 
are simple exclusion (SEP), zero range (ZRP) or misanthropes (MP) processes (see [1], [8]). 
Actually, these examples share two features: first, their infinitesimal transitions consist in 
the jump of a single particle (from a site x to another site y); they are also of misanthrope 
type, in the sense that their transition rates depend at most on x, y, and on occupation 
numbers at x and y. These limitations on the possible variations of the conserved quantity 
are at the root of the construction of basic coupling, through which attractiveness is shown 
to hold. 
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In this paper, we consider conservative particle systems of misanthrope type but allowing 
more than one particle to jump at a time from a given site to another one. This seemingly 
slight generalization requires to restart almost from scratch to check attractiveness through 
a coupling construction. We obtain necessary and sufficient conditions on the rates for at- 
tractiveness, giving an explicit construction of an increasing coupling (that is, which shows 
that the process is attractive). Using these results, we determine the extremal invariant 
and translation invariant probability measures of the process. Finally, we illustrate our 
approach on a few examples for which we derive hydrodynamics. 

The conservative particle systems we study are Markov processes {rjt)t>o defined on a 
state space VL = X^, with X a subset of Z, and S = Z"^; rjt is the configuration at time t 
and f]t{z) its value at site z ^ S. Depending on X, one of the two following interpretations 
can be used: for X C N, rit{z) is the number of particles at site z; ii X C Z, |?7t(-2)| is the 
number of unit charges at site z (either positive for rjtiz) > or negative for rit{z) < 0). The 
evolution of the process can be informally described as follows. Given a pair {x, y) G 5*^ 
and an integer A; > 0, with a rate ^^(x)ri{y)(y " ^ particles (respectively k positive 
charges) 'attempt a jump' from site x to site y. By this, we mean that the values at sites x 
and y are changed to r]{x) — k and r]{y) + A; as long as the resulting configuration belongs 
to the state space. The conserved quantity is ri{x) + ri{y), the total number of particles 
(resp. the total charge) involved in the transition. 

The SEP-ZRP-MP are examples of such processes, for which X C N, A; = 1 is the 
only allowed move and the rates take the form Tl^^i^y — x) = p{y — x)b{a,P) for a,/? G 
X, x,y & S, p{-) a translation invariant probability kernel on 5*, and b{-, ■) a function of the 
configuration values at the two involved sites. These processes (see [18], [1], [8] for details) 
are attractive, that is, the coordinate-wise (partial) order on configurations is preserved 
(stochastically) by the dynamics (see Section [2] for precise definitions), if and only if 

(A) The rate 6(a,/5) is a function of a, (3 non- decreasing in a and non-increasing in (3. 

Attractiveness for interacting particle systems is derived by constructing an increasing 
Markovian coupling. A coupled Markov process {C,t, Ct)t>o is a process on i7 x 1] such that 
each marginal is a copy of the original process {r]t)t>o- The coupling is increasing if partial 
order between marginal configurations is preserved under the coupled dynamics. 

The simplest candidate for this purpose is the basic (or Vasershtein) coupling. It is 
designed to make the two coupled processes {^t)t>o, (Ct)t>o 'as similar as possible', and its 
construction can be obtained by relying on three rules: 

(i) coupled transitions have the same departure and arrival sites for the two marginals, 

(ii) an identical change in both processes is given the largest possible rate, 

(iii) non coupled transitions are supplemented to get the correct left and right marginals. 

In the context of SEP-ZRP-MP, these rules give: by (i), the same probability p{y — x) 
is taken for coupled jumps from site x to site y. Then, if the values at sites x and y are 
^t-(x) = a, ^t-{y) = P, Ct-{x) = 7 and Ct-{y) = three changes may occur (remember 
that at most one particle can jump in a transition): 

* with rate (6(a, P) A 6(7, S))p{y — x), a particle jumps from x to y simultaneously in both 
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processes; 

* with rate (b{a, [3) — b{a, [3) A 6(7, S))p{y — x), only a .^-particle jumps from x to y; 

* with rate (^(7, 5) — b{a, (3) A 6(7, S))p{y — x), only a ^-particle jumps from x to y. 

By rule (ii), the last two changes are mutually exclusive, and their rates are fixed by 
rule (iii). Notice that the resulting Markovian coupling is not necessarily increasing. In 
fact, there are exactly two instances in which the coordinate- wise order ^t- < Q- might 
be destroyed, namely (a < 7 and P — S) or (a — ^ and P < S). In the first (resp. second) 
one, the isolated jump of a ^-particle (resp. a ^-particle) would break the order between 
configurations. Attractiveness requires that the corresponding rates are equal to zero, or 
equivalently 

b{a, 6) < 6(7, 6) for all 5 and a < 7 
fe(7, l3) < 6(7, S) for all 7 and /? > (5 

which is exactly the content of Condition (A). 

Beyond attractiveness, the basic coupling construction turns out to be essential to char- 
acterize the set (Tn<S)g of extremal invariant and translation invariant probability mea- 
sures of {r]t)t>Q, through a control of the evolution of discrepancies between the marginals 
of (^t,Ct)t>o- There is a discrepancy at ,x G iS at time t if ^t(x) ^ Ct(^)- SEP-ZRP-MP, 
in a coupled transition the number of discrepancies on the involved sites remains con- 
stant whenever the values of the two marginal configurations are ordered, but decreases 
otherwise. This is a consequence of the fact that Condition (A) (applied twice) leads to 
inequalities between rates on unordered configurations, 

b{a, P) < 6(7, S) for all q; < 7 and P > S 

This decreasing property is crucial: under suitable conditions on p{y—x) and b{a, P) , it per- 
mits the identification of (Xfl S)^ as a one parameter family {^p}p of product probability 

measures, where the parameter p fixes the average particle density per site. 

Beyond these classical examples, the requirements needed for attractiveness through 
rules (i-iii) might be so drastic that the basic coupling construction is not possible. We 
therefore have to weaken those rules in order to get attractiveness through a coupling 
construction imposing less restrictive conditions on the rates. At the same time, we have 
to ask for a decrease of discrepancies, which came for free within basic coupling but could 
be less obvious in a larger setting. 

Our starting point will be the more investigated case of pure jump processes with denu- 
merable state space, which provides us with a set of necessary conditions on the rates for 
attractiveness. Restricting ourselves to the class of conservative systems we consider in this 
paper, it turns out that it is possible to construct a Markovian coupling, which is increas- 
ing whenever these necessary conditions are fulfilled, hence proving also their sufficiency. 
Another property of this coupling is that discrepancies are non increasing. However, a de- 
crease of discrepancies requires additional irreducibility conditions. Under them, we prove 
that (Xn iS)g is a one parameter family of probability measures. This gives the first step 
in the derivation of hydrodynamic limits under Euler scaling, which we explicitely achieve 
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for two one-dimensional examples, where an alternative proof of macroscopic stability is 
needed in order to follow the approach of [1]. 

It is important to notice that the restriction to single particle jumps {k = 1) is essential 
in the basic coupling construction and in the derivation of attractiveness, in particular 
in the two following instances: it induces implicitely a total order on 5* and a partial 
lexicographic order on the state space; it drastically restricts the ways by which the order 
between configurations can be broken. These two points make the weakening of basic 
coupling rules useless and hide the interplay between partial order, conservation rule and 
optimization of coupling rates, though its understanding seems necessary to deal with more 
complex systems. Working with k > 1 gives some hints in these directions. In a forthcoming 
paper, we will consider attractiveness for conservative particle systems with speed change, 
a problem which requires a fine tuning of these (yet mostly hidden) properties. 

In the last part of this paper, we will develop our results on various examples: First, 
in the well known cases where only k = 1 is possible (thus in particular SEP-ZRP-MP) 
our construction reduces to basic coupling. Moreover, our irreducibility conditions relax 
restrictions on the rates for ZRP-MP and allow to determine of (X fl iS)^ in more general 
cases than previously treated. 

As soon as we depart from k < 1 case, both constructions show very different features: we 
focus on two one- dimensional examples, a discrete nearest neighbor stick process in corre- 
spondence with a generalized (that is, with nearest neighbor jumps) discrete Hammersley- 
Aldous-Diaconis (HAD) model, for which k is unbounded, and a two-species asymmetric 
exclusion process with charge conservation, for which k < 2. A continuous stick process 
with totally asymmetric jumps was introduced in [20] in relation with the (continuous) 
HAD model, and the discrete HAD process was studied in particular in [TUj. The two- 
species asymmetric exclusion process was introduced in [H] in the context of solid-on-solid 
interface dynamics and independently in [22j, [23]. Particular versions of this dynamics 
have been investigated in [12], [TTj in order to derive hydrodynamics in non-attractive 
cases. 

The paper is organised as follows. In Section[2l we first recall the definition and properties 
of attractiveness. Then we introduce the conservative models {rjt)t>o we are interested in, 
and state our main results: the set of necessary and sufficient conditions on the transition 
rates of {rit)t>o for attractiveness; the construction of the Markovian increasing coupling it 
involves, and under which discrepancies decrease. Proofs are given in Section [31 In Section 
we obtain (TDS)^ for (?7t)t>o under additional irreducibility conditions. In Section [71 
we analyze the hydrodynamic behavior of {'rit)t>o when the dynamics is one-dimensional 
and nearest neighbor. We have illustrated our results on examples, the classical ones, SEP- 
ZRP-MP, as well as on the stick process and the two-species asymmetric exclusion process: 
In Section [H and Appendix (Section [8]), the coupling construction is presented in some 
detail and in Section [6] we determine (X fl S)^ for each example. In section [7] we obtain the 
hydrodynamic limit of the stick process and the two species asymmetric exclusion; for the 
latter we prove heuristic results given in [9] and rederived more recently by [22], [23] . 
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2. Preliminaries and results 

2.1. Background. We first recall the general definitions and theorems about attractive- 
ness and couplings that we use in this paper. 

Let {r]t)t>o be an interacting particle system of state space fl = X'^, with X a subset 
of Z, and S = Z'^. We denote by T(t) the semi-group of this Markov process and by C 
its infinitesimal generator. We mainly refer to [TS], which deals with Feller processes on 
compact state spaces. Since we also consider other cases, we assume that 

{Vt)t>o is a well defined Markov process on a subset fig C ^2 
/r, -|\ such that for any bounded local function / on Qq, 

v„.n„.lmjm^ = £/(„)<+oo 

A coupled process {^t, Ct)t>o of two copies {^t)t>o and {(t)t>o of {rjt)t>o is a Markov process 
with state space VLq x VLq, such that each marginal is a copy of the original process {rjt)t>Q- 

There is a partial order (coordinate-wise) on the state space f2, defined by 
(2.2) Ve,CG^^, e<C^(Va;G5, ^(x) < C(a;)) 

Definition 2.3. (see [15]) Let W he a set endowed with a partial order relation. A subset 
V cW is increasing if 

(2.4) V/ G V, m e ly, I <m^meV. 
A subset V cW is decreasing if 

(2.5) \fl eV,meW, l>m^meV. 
A function f -.W -^M. is monotone if 

(2.6) V/, m G ly, l<m^ f{l) < f{m). 

Examples 2.7. (114] ) Let I G W . Then Ii = {m & W : I < m} is an increasing set, and 
Di = {m & W : I > m} a decreasing one. 

Remark 2.8. For any subset V C W, 

V is increasing <^=^ W \ V is decreasing <^==^ ly is monotone. 

We denote by M. the set of all bounded, monotone continuous functions on Q. The 
partial order (12. 2p induces a stochastic order (see [15]) on the set V of probability measures 
on Q endowed with the weak topology: 

(2.9) Vz/, u'eV, u<u' ^{yf EM, uif) < u'if)). 

Remark 2.10. v < v' <^=^ ^iY) ^ ^'(Y) for all increasing subsets V dVt. 

Indeed, on the one hand use Remark WT^ and on the other hand, notice that for any non- 
negative f E Ai and all s > 0, the set {/ > s} is increasing, and z/(/) = z/{/ > s}ds. 
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Theorem 2.11. flEl Theorem II. 2. 2] For the particle system {rit)t>o, the following two 
statements are equivalent: 

(a) f eM implies T{t)f G M for all t > 0. 

(b) For eV, u < u' implies vT{t) < v'T{t) for all t > 0. 

Definition 2.12. pLSi Definition II. 2. 3] The particle system {r]t)t>o is attractive if the 
equivalent conditions of Theorem \2.11\ are satisfied. 

By Remarks 12.81 and I2.10[ it is enough to check Theorem I2.1ir a) for indicator functions 
of all increasing sets. 

In the following our main object is the infinitesimal generator C of the process since we 
look for necessary and sufficient conditions on the transition rates that yield attractiveness 
of the particle system. As mentioned in [15], this was investigated by many authors for 
pure jump processes (see also [13], [211 Chapter 4]) on a countable state space. In such 
cases, sufficiency consists in writing the semigroup as an exponentiation of the generator 
to check Theorem 12.11( b). This follows from Strassen's Theorem, which links stochastic 
order to coupling (for a thorough analysis of this link, see |15j): 

Theorem 2.13. [HI Theorem II. 2. 4] For eV, a necessary and sufficient condition 
for V < v' is the existence of a probability measure fi on Q x Q, called a coupled probability 
measure, which satisfies 

(a) /i{(^, Q:^eA} = z/(A), /i{(^, Q : ( e A} = u\A), for all Borel sets AeVt, and 

r^;M(e,c):e<c} = i. 

Here we take the usual but longer route to derive attractiveness for interacting particle 
systems (see [IB]), and check Theorem 12.11( a) by constructing a Markovian increasing 
coupling, that is a coupled process (^j, Ct)t>o with the property that ^ Co implies 

(2.14) p(€o,Co)|^^ < C*} = 1 

for all t > 0, where P^^o,(o) denotes the distribution of i^t,Ct)t>o with initial state (^oCo)- 
Our strategy for proving and using attractiveness will consist in three steps: 

(a) obtain necessary conditions on the transition rates of the model; 

(b) construct a Markovian coupling and show that it is increasing under conditions obtained 
in (a), thus proving that they are also sufficient; 

(c) verify that discrepancies cannot increase. When supplemented with some conditions 
on the rates, this will allow to determine (Xfl S)^. 

2.2. The model, attractiveness and coupling. From now on, we restrict our analysis 
to the following class of models. 

Let S = Tj"^ he the set of sites and X C Z be the set of admissible values on each site. 
In our examples, X will be either a finite subset of Z or N. The infinitesimal generator C 
of the process {rjt)t>Q onVt = is given, for a local function /, by 

(2.15) cf{r^) = J2 E x::;{v)Y.^Uy-^)(f(siyv)~f{v)) 

x,y&Sa,P&X feeN 
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where Xx',y indicator function of configurations with values on {x,y), 

(2.16) x'iM = {l .f';W="and,fe) = ;3 

I U otherwise 

5*^^ is a local operator performing the transformation whenever possible (the value k = 
is not excluded) 

{r]{x) — k li z = X and ri{x) — k G X, r](y) + k E X 
Viy) + k if z = y and vix) -ke X,r,{y) + keX 
T]{z) otherwise 

This particle system is conservative, with ri{x) + ri{y) the conserved quantity in a transition 
between sites x and y. 

Unless X is finite and the rates F' (.)'s have finite range (in which case we refer to [T8l 
Chapter 1]), additional conditions on the transition rates, and/or a reduction of the state 
space, are required to ensure that fl2.15p is the infinitesimal generator of a well defined 
Markov process. Since such conditions differ for the specific model one deals with, we 
assume (12. ip and state here only a common necessary restriction on the rates (complete 
precise assumptions being given on examples). For all z E S, a, P E X , the rates T'^^z) 
satisfy 

(2.18) J^T^Jz) <oo 



a,l3\ 

ken 



For notational convenience, we will often drop the explicit dependence on z of F^^^(z) ; 
we also set 

(2.19) T'l^p = 0iia-k^X,(3 + k^X, or = 0. 

Definition 2.20. The notation {a, (3) < (7,5) is equivalent to a < '-f, /3 < 5; the two 

pairs (a,/5),(7,5) are ordered if (a,/3) < (7,5) or {a,j3) > (7,5); they are not ordered 
otherwise, that is when (a < 7,/3 > 5) or (a > 7, /5 < 5). 

The main result of this section is 

Theorem 2.21. The particle system {rjt)t>Q is attractive if and only if for all {a,P,'j,6) G 
X^ with {a,/3) < (7,5), for all {x,y) G , 

(2.22) V/>0, J2 < T.<s{y-x) 

k'>5-f3+l l'>l 

(2.23) VA;>0, ^ r^;^(^ _ a;) > F!;,,(y - x) 

k'>k Vy-f-a+k 

Note that by (12.181) the above sums are finite. We now give the succession of steps which 
leads to this theorem. Proofs are postponed to Section [3J 

First, we prove necessary conditions for attractiveness: 
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Proposition 2.24. // the particle system {i]t)t>o is attractive, then for all (a, (3, 7, 5) G 
with {a, (3) < (7,5), for all {x,y) e S^, inequalities f lZ^ - fl^^ hold. 

The other part of Theorem 12.211 is harder and is obtained by constructing explicitly a 
Markovian coupling, which appears to be increasing under conditions (I2.22p - fl2.23p . The 
evolution of the coupled process (^j, Ct)t>o we look for is defined through its infinitesimal 
generator £, whose rates will be derived from those of C Here, as in the basic coupling 
construction, the same departure and arrival sites x and y are chosen for jumps of coupled 
particles. 

Proposition 2.25. The operator C defined on VLq x VLq as 

^/(^'O = E E E ^"f X 

a;,yG5 o,/3gX 7,(5gX 

(2-26) E G^!pn^y - ^) (/(^' ^U) - /(^' 0) 

k,l 

with coupling rates given by (12.271) - fl2.30p below for all (a, /5, 7, 6) G X'^ and all non- 

negative k, I as functions of the initial rates F^'^ and Fjjl^, is the generator of a Markovian 
coupling between two copies of the Markov process defined by f l2.15p . 
For all positive k, I, we set 

(2.27) A(r!^,, - r;,, A (s^,^ - s^,^ A j:\^s)) 

(2.28) Gj;,^^,, = t\^,-t\^,a{j:%-j:%ai:\^s) 

(2.29) = r'^.,,-ripA{Ks-Kp^Ks) 

(2.30) G'^X.^s = 

where the S^^ 's denote partial sums of rates for k > 0, 
(2-31) K,, = E C- 

k'>k 

In other words, at time t, configurations C,t-, (t- are changed into = S^^^y^t-, (t = 
Sl yCt' with rate G'^'[^^-^ ^ -{yyc -(x) c - {y)^y ~ -^^^^^ with superscript '/c; 0' (resp. '0; /') 
correspond to uncoupled changes in the ^. (resp. C) process. 

Since the coupling rates are constructed for a fixed pair of sites {x,y), we dropped from 
the notation the explicit dependence of the rates on {x,y) and wrote F^^^ for T^j^{y — x) 

and ^51^.^ ,5 for GjJ^.^^^d/ - x). 

Remark 2.32. The partial sums S^^'s in (I2.3ip are well defined by (12.180 . As a direct 
consequence of expressions (12.270 - (I2.30p for the coupling rates, the coupled process is well 
defined on ^Iq x (cf. (12. ip ). 
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Remark 2.33. Formulas (12.271) -f l2.30p are symmetric under exchange of marginals: for 
all (a, P, 7, 5), k and I 

/^k\ I k 

Remark 2.34. Like in the basic coupling construction, the diagonal coupling rates are 
equal to the original rates of the marginal processes: for (a, /?, k) = (7, 6, 1), (12.27^ becomes 



/^k;k / pfc T^k A ^'V^ A V'^ ^ I P 



k 

a,0 



In the next two propositions, we give some hints on the construction and structure of 
the couphng rates; for any fixed (a, /?, 7, 6) G X^, non zero couphng rates are located on a 
"staircase-shaped" path in the {k, I) quadrant: most of the couphng rates are zero and for 
n > there is at most one pair (/c, I) such that k + 1 = n and G^'^.^^ > 0. Ahernatively, 
the couphng rates can be defined by induction. 

Proposition 2.35. For any 7,5) e X'^, all z e S, let (/co,/o) = (0,0) and Va,i3;'y,5 = 
{{ki, li)}i>o be the sequence of points inN x N defined for all i > as 



(2.36) 



ki — ki^i + 1, li — h-i if — ^7,5 

ki = /cj-i, k = U-i + 1 otherwise 



The coupling rates defined by equations (12.271) - (12. 3Up satisfy 
(i) For all (kiJi) G Va,f3;-t,s 

(2.37) Gjl^;^,, = sj- V si;- - sj^ v j:^, 

(a) For all {k,l) ^ 'Pa,f5-'y,& 

(2.38) G'^!pn,s = ^ 

Proposition 2.39. The set of coupling rates defined by equations (I2.27P -( 12.301) is the 

unique solution of the recursion relations: 

For all (a, /9, 7, 5) G X^, all z & S and all positive k,l, 

(2-40) Gji,,, = K,-Y.G'Sn,s)^K^-Y.<Ls) 

l'>l k'>k 

(2-41) ^a,/3;7,<5 = ^\.S ~ ^.5 

k'>0 



(2-42) = Ti,-Y.G'S 



7,5 

Z'>0 



(2.43) G'^' = 



a, (3; 7,(5 

We end this subsection with a proposition giving the last step to Theorem 12.211 
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Proposition 2.44. Under conditions fl2.22p - fl2.23l) . the Markovian coupling defined in 
Proposition \2.25\ is increasing. More precisely, if {a, (3), (7, 6) are ordered, then the coupling 
rate G^^.^g (with non-negative k,l) is nonzero only if the increments satisfy 

.2 U-ke {-{5 -/?),■■■, 7 - «} ^fia,/3)< (7, 6) 

^' ' \l-kE{-ia-^),--- ,P-S} z/ (a, /3) > (7, 5) 

2.3. Evolution of discrepancies. Attractiveness expresses that two ordered configura- 
tions remain such (stochastically) under the dynamics. Many applications require also a 
control over unordered pairs of configurations. This can either be a consequence of attrac- 
tiveness or require some additional hypotheses, that we investigate in this subsection. 

Definition 2.46. In the coupled process {C,t, Ct)t>o, there is a discrepancy at site z E S 
at time t if ^t{,z) 7^ C,t{z). This discrepancy is positive (resp. negative^ of width a > if 
^t{z) - Ctiz) = a (resp. ^tiz) - Ct{z) = -a). 

We first show that under the Markovian increasing coupling constructed in Proposition 
12. 25^ the sum of the widths of the discrepancies (involved in a transition) between the two 
coupled processes {^t)t>o and (Ct)t>o is non-increasing (later on, we will sometimes forget 
the word 'width', and speak of non-increasing or decreasing discrepancies). 

We fix X, y two sites of 5*. Denote by Cx,y the generator of jumps between them, 

(2.47) W(e,C)= E E ^"f (Ox2;S(C) X 

k,l 

+ E - y) (/(^Le, 5j,.C) - /(e, C))) 

k,l 

for a local function /, and two configurations ^, C ^ ^o- 

Lst f^y = fy :^ be the functions which measure the width of positive or negative discrep- 
ancies between ^ and ( on sites x and y 

(2.48) /,±,(e,c) := ^.^(e,c)+^j(e,c) 

The quantities 

(2.49) A^{ax),ay)X{x)X{y),kj) := f^is^J^ S^yO ' f^,yi^X) 

measure how, in a coupled transition on sites x, y, discrepancies evolve. Since by the 
conservation rule, 

C) - f7,yi^, = ia^) + ay)) - (C(^) + Ciy)) = f^^'j^ ^1/) - f7MJ^ ^i/) 

the two quantities A^*" in fl2.49p are equal and we set 

(2.50) A{ax),ay)X{x),ay),k,i) = A^{ax),ay)X{x)X{y),k,i) 
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with the meaning that positive and negative discrepancies disappear (or get reduced) by 
pair annihilation. Notice the symmetry 

(2.51) A(e(y), e(x), ciy), cix), I, k) = A(e(x), e(z/), ax), ay), k, i) 

because by exchanging respectively ^(x) and ^(y), ({x) and ({y), k and /, 

my) - 1) - (ay) - k)]^ = my) + k) - (ay) + 1)]^ 

mx)+i)-{ax) + k)]^ = mx)-k)-{ax)-i)v 

Next theorem details how discrepancies decrease. 
Theorem 2.52. For all (^, C) e Qq x ^o, {x,y) e S'^ 

(2.53) :c.,./j:,(e, C) = c.j-j^, C) < 

More precisely, setting (a,/?, 7, 5) = (^(2:), C(?/), C(a^), C(y)); the coupling rate G^^^p.^^^{y-x) 
may he positive only for non negative increments k, I satisfying: 

(2.54) - [a - 7]+ - [5 -/?]+</- A; < [7 - a]+ + [/3 - 5] + 
and in that case 

A(a, /3, 7, 5, A;, /) = for (a, /?), (7, 5) ordered 
A(a;, /3, 7, 5, /c, /) = /or (a, /?), (7, 5) not ordered 



(2.55) 



andk-l e {0, (a - 7) + (5 - /?)} 
^ A(a, 7, 5, fc, /) < otherwise 



In ( 12.551) . the first line expresses attractiveness. The second line deals with discrepancies 
of opposite signs and contains two cases: when k — I = 0, the same number of particles 
moves together in both processes so that discrepancies do not change (as in basic coupling, 
for which k = I = 1); when k — I = (a — 7) + (5 — /5), the positions of discrepancies at 
sites X and y are exchanged. The third line expresses a decrease of discrepancies, it will 
be the key tool in the determination of {Xf\S)^ (see Section [5]). Notice that for basic 
coupling, this case corresponds only to uncoupled jumps, while here it may involve non- 
zero values for both k and /; the jump gets the two marginals closer by merging together 
two discrepancies of opposite signs. The resulting pairs of values need not be ordered since 
discrepancies can reduce their widths without disappearing and possibly exchange their 
signs (as in the second line). This new phenomenon of (partial) exchange of discrepancies 
does exist beyond basic coupling (it requires \k — l\ > 2) and induces some changes in the 
behavior of the coupled process. More precisely 

Definition 2.56. The increasing coupling defined by (12.261) allows exchanges of discrepan- 
cies whenever there exist {x,y) G S"^, (a,/3,7,(5) G with (a — 7)(/3 — 5) < 0, non-negative 
k,l with \k — l\ > |a— 7|V|/5— 5| such that C'^p.^ ^{y—x) > 0. The exchange of discrepancies 
is total if k — I = (a — 7) + (5 — /3) and partial otherwise. 
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Lemma 2.57. The increasing coupling defined by fl2.26p does not allow exchanges of dis- 
crepancies if and only if for any {x,y) G S'^, (a,/5,7,5) G X^, non-negative k,l such that 

^a!f3;'y,6iy - ^) > ^ > ^C haVC 

(2.58) - {[a - 7]+ V [5 - /3]+) < / - A; < [7 - a]+ V [/3 - 5] + 

Remark 2.59. When (a, /5), (7, 5) are ordered, (12.581) reduces to inequalities (12.451) for 
attractiveness. Otherwise (12.581) is stronger than (I2.54p (which is a consequence of attrac- 
tiveness). 

A first application of Tlieorem 12.521 is to determine extremal invariant and translation 
invariant probability measures for the model (see Theorem 15.131) under irreducibility con- 
ditions on the rates, that cope with exchanges of discrepancies and other phenomena which 
appear beyond the case k < 1. 

A second application is macroscopic stability, an essential step to derive the hydrody- 
namic limit of one-dimensional models (see Section [7]). This property, introduced in [71 
Proposition 3.1] for one-dimensional finite-range simple exclusion processes, gives a con- 
trol of the time evolution of discrepancies in the coupled process. However, its proof does 
not extend to models with jumps of size > 1, in particular because the latter induce 
exchanges of discrepancies. So we prove here 

Proposition 2.60. We assume that the particle system defined by (I2.15P is one dimen- 
sional (S = TL) with only nearest neighbor transitions, and that the increasing coupling 
defined by (12.261) does not allow exchanges of discrepancies. Then, if the coupled process 
i^s, Cs)s>o is such that '^^^zl^oi^) + (oi^)] < +00, we have for every t > 0, 

(2.61) 5(6,0) < 5(^0, Co) 
where, for x G Z, 

(2.62) s^i^sXs) ■.= J2iUy)-Cs{y)], Si^sXs) := sup |s,(es,C)| 



Remark 2.63. Both restriction to nearest neighbor interactions and Condition ([21 



are necessary to get the macroscopic stability property (12.611) . See Section\^for counter- 
examples. 

Proposition 2.64. The following inequalities imply (I2.58p .- For all {a,P,j,6) G X"^, non- 
negative k, I, 

(2.65) J:\^s > L = /+[a-7]+V[5-/3] + 

(2.66) T^Ip > K = k + [^-a]+y[/3-6]+ 

Remark 2.67. When (a,/5),(7,5) are ordered, (IZeSj) - (12361) coincide with (12:221) - fl2:23|) 

and thus reduce to attractiveness conditions. 
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3. Proofs. 

Proof of Proposition I2.24L 

We proceed in two steps. The first one yields an inequality for the generator of an attractive 
system, when apphed on indicators of increasing sets. Such inequahties were first derived 
in [19] for Markov jump processes on a countable state space E, with an infinitesimal 
generator bounded in li{E). The second step specializes the previous inequality to our 
model on which it reads (^^ - (K23\f . 

Step 1. Let C) G 1^0 X ^0 be two configurations such that ^ < C- Let V dVthe sji 
increasing cylinder set. \i &V oi C, , 

(3.1) V(0 = ly(C) 
By attractiveness, for all t > 0, 

™v)(o<mt)M(c) 

(use Remark 12. 101 and Theorem 12.11( a)). Combining this with (13.11) . 

t''mt)iym - v(0] < t-\{T{t)iym - v(c)], 

which gives, by Assumption (12.11) . 

(3.2) (£V)(0 < (£V)(C) 
We have 

(^M(0 = E ^tv^^^ll^'^dy - ^) {MsU) - MO) 

x,yeS a,l3eX k 

= -v(oE E ^",f(^)Ew(^i.orUz/-x) 

x,y£S a,l3€X k 

(3.3) +ino\y(o E E xziioY.'^vKyO^lAy - ^) 

x,yGS a,f3eX k 

Step 2. We fix {x,y) G 5^, x ^ y, (a,/?, 7,5) G X"^, with (a,/3) < (7,5), and two 
configurations (^,C) in x such that C,{x) = a, C,{y) = (3, C,{x) = 7, ({y) = 6, and 
^(z) = ^(2;) for all z ^ x,y. Thus ^ < C. 

We fix numbers {pz G X,z G 5"), which satisfy Px < C,{x) = a, Py > ({y) = 6, and 
Pz = ^z) = ({z) for all z G 5 \ {x,y}. We denote by || ■ || the Li-norm on 5*. For 
n > \\y — x\\, we define 

Iy{n) = {t] E Q : ri{z) > Pz, for all z 7^ x such that \\y — z\\ < n} 
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By Example \2.7\ Iy{n) is an increasing cylinder set. Since neither C, nor belong to Iy{n), 
we have using fl3.3l) 

x',y'(^S a',l3'eX k' 

fc'GN z:\\y-z\\>n fc'eN 

and similarly 

I'eN z:\\y-z\\>n /'SN 

SO that (13.21) writes 

E E E 

i{/3+fc'>py}r^(2)^^(y — 2:) 

k'€N:f3+k'>py z:\\y-z\\>n fc'SN 

< E Ks(y-^)+ E E i{^+''>Mrc(.),5(i/ - ^) 

Z'GN:<5+/'>py 2:||s/-z||>n /'GN 

Taking the monotone limit n —>■ 00 gives 

(3.4) E r^>-^)< E <s{y-x) 

k'eN:f3+k'>py I'eN: 5+l'>py 

A similar argument can be used with the decreasing set 

Dx{n) = {t] E Q : r]{z) < p^, for all z ^ y such that ||x — z\\ < n} 

Configurations ^ and C belong to its complement QQ\Drc{n) (which is increasing by Remark 
12. 8p and inequality (13. 2p for this set leads to 

(3.5) E <(siy-^)> E <s{y-x) 

k'eN:a-k'<px I' GNi-y-l' <px 

Finally, taking = 5 + / + 1 in ([S3D and = a - A; - 1 in ([S3]) gives ^Mj-^M)- ■ 

In order to prove Proposition 12.251 we give equivalent expressions for the coupling rates. 

Lemma 3.6. The following expressions are equivalent to (12.271) -( 12. 29p for all [a, (3, 7, 6) G 
and positive k, I: 

(3-7) Gj;^^,, = r^,^ A - A E^J) - r^,, A (S!^,, - E^,, A E^,,) 

(3-8) GJi^,, = (E^7/ - E^7,i A e!^,,) a r;,, - (E^,, - E^,, A e!^,,) a T\^, 

(3-9) Gj;,^,,. = KJ - a KJ) - - Kf^ A e;,,) 

(3.10) G^^,,, = {E'J - A E^,,) - (E^,, - E^,, A E?,,) 
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Proof. We prove only fl3.8p . (13.101) : the two other equations follow by symmetry from 
Remark [2.331 Let (a,/5,7,5) G X"^, k> and / > 0. Using the elementary identity 

(3.11) Va, 6,c>0, aA{c-bAc) = {a + b)Ac-bAc 

with a = ri^^, b = T^lp, c = leads to (we have T^^ = l^J - S^^^) 

- A (S;,, - a S;,,) = [{a + b)-{a + b)Ac]-[b-bAc] 

(3-12) = [KJ - ^'J A S;,,) - (S^,, - S^,^ A S!^,,) 

Recalling (12:291) and setting / = in (l3J2ll gives (l3A0ll . Inserting (l3A2l) in (E^ZD, and using 
identity (EH]) again with a =(E^-i-E^-iASi^,,)-(S^_^-S^,^AEi^,,), 6 = E^^^-E^^^AS^^^, 
and c = Ti^^^ gives 

a[T\^s'K,s^Kp-Kp^Ks)] 

= i^'J - ^'J A e!^,,) a r!^,,- (E^,^ - E^,^ A e!^,,) a r!^,, 

which is (I33D. ■ 
Proof of Proposition I2.25L 

Let (a, /?, 7, 6) G and positive k, I. Using (12.291) . (13.71) and a telescopic argument, partial 
sums of coupling rates read 

I I 

E^k;l' _ ^k;0 ^k;l' _ -pk A ('v' V'^ AV' ^ 

^a,(3;'y,S — '^a,l3;'y,S ~^ Z_^^a,P;'t,5 — ^ a,l3 ^ a,l3 \^'y,5 ^a,f3 ^■y,5) 



l'=0 l'=l 



In the limit / oo, E^^ ^^ — > by Assumption (12.181) . thus 
(3.13) Inn r^,^ A (E^,, - E^,^ A E^,,) = 

This gives the correct left marginal for all A; > 0, 

-\ A\ \ ^ ryk;l' ^k;0 , \ ^ ^k;l' -pA 



^k 

a,l3 

Z'>0 l'>0 



The right marginal can be treated in a similar way using (12.281) and (13.81) . We get 
(3.15) hm (E^,^ - E^,^ A E^ A r;,, = 



fc^oo 

thus, for alH > 0, 



(3.16) ^ ^1/3; 7,5 - ^a,^;7,5 + ^t,'l3;^,S " ^7, 



7,5 

fe'>0 k'>0 



Proofs of Propositions 12.351 and 12.391 rely on the following lemmas. 
Lemma 3.17. The coupling rates defined by (I2.27P - (12.301) satisfy the recursion relations 
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Proof. Equations (I2.4ip - fl2.42p are verified by fl2.27l) - fl2.29l) as a direct consequence of 
dSHD-dSlS]) (and fl^:^ is identical to ([233D)- Equation fl^iU]) requires a partial resum- 
mation of the coupling rates : using (13.71) . (13.131) on the one hand and (13. 8p . (I3.15P on the 
other hand gives 

(3-18) E <Ls = r',, A (s;,, - s^,, A s;,,) 

l'>l 

for > and / > 0, and 

(3-19) E GS^,,, = (Si, - A S^,) A 

k'>k 

for /c > and / > 0. Both right hand sides above appear in (I2.27P : substituting them by 
the respective left hand side of fl3J8l) and fl3J9D in fOTD gives ^M). ■ 

Lemma 3.20. For all (a,/3,7,5) G , any solution of (12. 40p - (12.430 satisfies for all 
{k,l) G N X N, 

(3.21) > S^, ^ VA:' < k, V/' > /, = 

(3.22) < S^, ^ VA;' > k, V/' < /, = 

Proof Inserting ([23T])-([232]) into fl2:i0|) writes 



0<l'<l 0<k'<k 



Substracting G^^.^^ on both sides gives 

(3-23) (E gS;.,)a($:gsl..)=o 

0<fc'<fc 0<i'</ 

On the other hand, we have 

k'>k]l'>0 k'>0]l'>l 

(3-24) = E gS;,,- E 

k'>k;l'<l k'<k;l'>l 

Let > S^^ ^. Suppose that the first sum in (I3.24p is zero; by non- negativity of the 
coupling rates, all terms in the second sum in (13.240 are zero, which gives (I3.2ip . Now 
suppose that the first sum in (13.241) is not zero. Then there are k* > k, I* < I such that 
Ga/3-7<5 > 0- Thus, for all I" > I, we have 



0<l'<l" 
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and by (I3.23p applied to each /" 

0<k'<k* 

By non-negativity of the couphng rates, all the elements in the above sums are zero; in 
particular, since k < k*, (13.21 1) holds. 

We skip a similar derivation for f l3.22p . I 

Lemma 3.25. The rates defined by (12.371) - (12. 38p are the unique solution of the recursion 
relations (ICTj) - (ICTD . 

Proof. It consists in an explicit resolution of ( ]2.40p -( l2.43p . We first prove that any solution 
is zero outside Va,i3;^,8, that is ( 12.380 . Let {k,l) ^ 'Pa,i3;'y,s and define i = k + I and 
{ki,li) as in ( 12.360 . Notice that by construction ki + U = i (hence k ^ ki, I ^ li), 
{ko, ki, - ■■ ,ki} = {0,1, - ■■ , ki} and {Iq, /i, ■■■,/»} = {0, 1, ■■■ , U}. 

Suppose first that k < k^. Let i' denote the greatest element of {j > : kj = k}. 
Hence = /cj/ + 1 and by (I2.36p . we have S^*'^ > s!^''^. Since k = ki' and I > k > h', 
GJi,,, = by (132ID. 

Suppose now k > ki. Then / < U and we denote by i" the greatest element of {j > : 
Ij = I}. Hence /j/'+i = k'l + 1 and by (12.360 . we have S^*'^ < e!^''^. Since k > ki > kin and 

I = h", Gj;,^^^, = by m^. 

We now are able to compute exphcitely the values of the rates on Va,f3--y,s- 

Take i >0 and {hji) as in ( 12.360 . Suppose first that ^^,0 — ^7,5- Then kj > ki for all 
j > i and thus, using (12.380 . 

j>i k>ki l>0 k>ki 

Suppose now that ^a,i3 < ^7,5- Then Ij > li for all j > i and we have 



^7,5 

j>i l>li k>0 l>li 

Putting both cases together gives for alH > 

j>i 

These relations lead to the expression (I2.37P : for all i > 0, 

j>i-l j>i 
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Proof of Propositions [2735] and [27391 

From Lemma 13.171 the coupling rates defined by fl2.27p - fl2.30p verify tlie recursion rela- 
tions fl2.40l) - fl2.43l) . which, by Lemma 13.251 have a unique solution given by 02.370 - 02.380 . 
Therefore expressions (^;^-I^Mj and l^3I^-I^M> are identical. ■ 

Proof of Proposition I2.44L 

Using conditions fl2.22p - fl2.23p . we now prove that the process is attractive. We need to 
show that = whenever < (7,5) and {a — k, P + k) ^ (7 — /, 5 + /). This 

last condition splits into two possible cases 

(3.26) l>-f-a + k 

(3.27) k>6-(3 + l 

In the first case, using fl3.26p . non-negativity of the F's and fl2.23p . we write 

l'>l l'>^-a+k k'>k 

Since S^^ ,^ < one has both ^ < ^a,p ^ind Sj^y/ < S^^^ and the two terms appearing 
between parentheses in expressions 03. 7p for k ^ or 03. 9p for A; = are equal to zero. 
Hence G'JJj;^,^ = 0. 

In the second case 03.271) . one writes in a similar way, using fl2.22p 

^a,/3 ^ ^a,/3 " ^ a,/3 - ^ - Z^ ^ 7,5 " ^7,5 

Inserting these inequalities into the expression for the coupling rate, that is 03. 8p for / > 
or 03.1 op for I = 0, one gets again G'^^.^ g = 0. 

In conclusion, when (a,/?) < (7,5), we have 

G^ai7,5 = 0if/-A:^{-(5-/3),---,7-«} 
The case (a,/3) > (7,5) is treated by symmetry using Remark [2.331 which yields 02.450 . I 
Proof of Theorem [2321 

We first prove that for all (a,/3, 7,5) G X'^, the coupling rate G^'^.^^ is zero if /c — / is 
outside some range depending only on a — 7 and P — 6. 

From Proposition 12.441 we already know that if (a, /5), (7, 5) are ordered, then fl2.45p 
is satisfied. We now show similar restrictions for non ordered pairs, for which equations 
f i2.22p - fi2.23p lead also, indirectly, to inequalities between partial sums of rates. Suppose 
for instance that a > 7, /5 < 5. We insert an intermediate pair of values and write 

(«,/?)> (7,/?) and (7,5) > (7, /5) 
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We apply equations (I2.22p . fl2.23p on both pairs of values and get two sets of inequalities 
valid for all > 0, 

k'>k V>k V>a-^+k k'>k 

fc'>5-/3+fc l'>k l'>k k'>k 

Combining them by pairs in order to eliminate the intermediate values (7,/?) gives 

Er?. < EC 

l'>k k'>k 
k'>5-l3+a--/+k l'>k 

which implies 

(3.28) K,s<K} ^ for alH > 

(3.29) T.l^^<T.lJ < Si^^forall A;>5-/3 + a-7 + / 

Inserting (13:281) in ([33D (or ([33]) if A; = 0) and (13:291) in ([3lD (or (l3:T0|) if / = 0), one gets 

(3.30) <^ai7,5 = if a > 7, /? < 5 and A; - / ^ {0, ■ ■ ■ , a - 7 + (5 - /?} 
We have also by symmetry through Remark 12.331 

(3.31) = if a < 7, /3 > 5 and / - A; ^ {0, ■ ■ ■ , 7 - a + /3 - 5} 

Thus, putting together (12.451) . (13.301) . (I3.3ip . the coupling rates G^'^.^^ are nonzero only 
if the increments satisfy 

'l-ke {-(5 - /5), ■ ■ ■ , (7 - «)} if a < 7 and /5 < (5 

l-ke {-{a - 7), ■ ■ ■ , (/5 - 5)} if a > 7 and /5 > 5 

I - k e {-{a - 7) - (5 - /3), ■ ■ ■ , 0} if a > 7 and P < 6 

J - k e {0, ■ ■ ■ , {'J - a) + {f3 - 6)} if a < 7 and (3 > 6 

All four cases can be written together as the single statement (12.540 or equivalently 

(3.33) |/ _ A; + l(a - 7) - \{[3 - S)\ < - ^\ + ^\(3 - 6\ 

These restrictions induce in turn a control on the sign of A(a, /3, 7, 5, /c, /) through the 
inequality: For all a,b,x in M such that |x| < \a\ + 

(3.34) [a + b + x]^ +[a + b-x]^ < 2[a]+ + 2[6] + 

Inequality (I3.34p is trivial if its left hand side is equal to 0. When exactly one of the two 
quantities on its left hand side is nonzero, (13.340 follows directly from the bound on x; 
otherwise, the value of the left hand side is independent of x and trivially bounded by the 
right hand side. 



(3.32) 
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We now apply (13.341) . choosing (in view of (13.331) ) 

(3.35) a = ^(«-7) b=^{[3-6) x = I - k + ^{a - ^) - ^{f3 - 6) 
we get 

(3.36) [{a -k)-{-f- /)]+ + [(/? + k)-{6 + /)]+ < [a - 7]+ + [/? - S] + 
Hence, recalling Definition (12.491) of A(a, /?, 7, 5, fc, /), and (12.511) . 

(3.37) G^;;,^ 7,5 7^ ^ ^(«' 7, 5, fc, /) = A(/3, a, 5, 7, l,k)<0 

To go from fl3371) to fl233|) . notice first that for A; = /, A(^(a;), ^(y), C(x), C(y), fc, = 0. 
We have, since /^(e, C) = f^^, (), and by (EUD, 

(3.38) J] (^'^7,5(1/ - ^) A(a, /5, 7, 5, fc, + G'J,t;5,7(a; - 2/) A(/5, a, 6, 7, A;) 

= E E (^)^2;^(o X 

(3.39) E(^'i7,5(^ - ^) + ^At;5,7(^ - ^))^(«' ^' 7' ^' 

By (I3.37p . each term of the sum in the right hand side of (13.381) is non-negative. Since 
by (I230D, :C,,,/+,(e,C) =^x,y/.";,(e,C), (ESSD is satisfied. 

To derive (12.551) . we first check that equality holds in (13.341) either for a + b = ±(|a| + \b\) 
or for X = ±(|a| + 16|). Recalling (13.351) . the first case corresponds to ordered pairs of values, 
(a,/3) < (7,5) or {a,P) > (7,5), for all values of k and /. The other one gives two new 
possibilities for non ordered pairs of values: either I — k = OoTl — k = — (a — 7) + (/9 — 5). 
Notice that in this last case, the values of the discrepancies at sites x and y are exchanged. 

■ 

Proof of Lemma 12.571 

Let {x,y) G S"^, (a,/?, 7,5) G X^, non-negative k,l such that G^'^^.^^d/ — x) > 0. Either 
(a,/9),(7,5) are ordered and (I2.58P reduces to the expression (I2.45P for attractiveness; 
or (a,/5),(7,5) are not ordered and (I2.58p is Definition 12.561 of absence of exchanges of 
discrepancies. I 

Proof of Theorem 12.601 

We fix (^0, Co) e X2 with T^yezMv) + Co{y)] < +00. We prove fimi) by checking that 

(3.40) CS{^t, Ct) < for alH > 
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We have 



x,ye'E,\x-y\=la,l3£X j,5eX k,l 

(3.41) (sup|s,(^J^6,^LCJ| -sup|.,(6,Ct)|) 
We note that 

(3.42) Wz^[xAy,xyy), s,{S',Jt, S^yCt) = ^.(6, Ct) 

because, for z < xAy, jumps involve only sites x, y and because, for z > xVy, the dynamics 
is conservative (cf. (12.501) ). Due to the restriction to nearest neighbor interactions, ||/— x| = 
1, so that Sz may change only for 2; = x A 

'sMju 5i,,C*) = s.(6, Ct)+l-kifz = x = y-l 
SziS'^Ju S^yCt) = s.(6, Ct)+k-liiz = y = x-l 



(3.43) 



and we have 

csi^tXt) = EE E ^";'+i(^*)^^'+i(c*)E(^S;7,^(i)+^A^^^ 

x&Z a,l3eX y,SeX k,l 

{\sx{^t, Ct)+l-k\y sup \s,{^t, Ct) I - 5(6, Ct)) 

^ EE E ^";'+i(^*)^^'+i(c*)E(^S;7,^(i)+^At^ 

(3.44) [\s^(^^Xt)+l-k\ -Si^tXt)^ 

Now we prove that if the coupling is increasing and does not allow exchanges of discrep- 
ancies, all terms in the last line of (I3.44p are zero. Let (Ct, Ct) be such that ^(x) = a, 
6(x + 1) = P, Ct{x) = 7, Ct(a; + 1) = 5. For all transition rates with ^^'^^.^,5(1) + 
Gf3a-S'yi~^) > 0' apply Lemma [2.571 Indeed, either G„'^.^^(l) > hence we have 
(I2.58p . or G''^''^. ^ ^{—1) > 0, which implies inequalities (12.581) again (by symmetry the latter 
stay the same when exchanging a with /3, 7 with 6, k with I). They have the equivalent 
formulation 

(3.45) A (7 - a) A (/? - (5) < / - A; < V (7 - a) V (/? - 5) 
Recalling (I2.62p we have 

(3.46) s^-i(6,Ct) = ^^(^.Ct) +7 - «, s^+i(6,Ct) = Sx(6,Ct) 
Thus, by adding Sx{C,t, Ct) to each term of (I3.45P we get 

Sxi^t, Ct) A s^_i(6, Ct) A Sx+i{^t, Ct) < Sxi^t, Ct)+l-k< s^(6, Ct) V s^_i(6, Ct) V s^+i(6, Ct) 
which implies 

(3.47) \s,,{^t,Ct)+l-k\ < |s,(6,Ct)| V|s,_i(6,Ct)| V|s,+i(6,Ct)| <5(6,Ct) 

and finally, inserting (13^^71) in (^ZM, we get (K^ . ■ 
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Proof of Remark 12.631 

We exhibit two counter-examples for which fl2.6ip is wrong: in (i), interactions are not 
only nearest neighbor; in (ii), fl2.58p is not fulfilled. For both cases, let (^o, Co) be such that 
for X < y, ^o{x) = a, ^o{y) = (3, Co{x) = 7, Coiv) = ^ with (a,/?), (7,5) not ordered, and 
^o(^) = Co(^) = for 2; ^ [x — l,y]. Let the first coupled transition occur at time t > 
between sites x, y at rate G'^p.^siv ~ 

(i) Assume that y = x + 2, ^q{x + 1) = 5, ^0(2; + 1) = 7, ^0(2; — 1) = Co{^ ~ 1) = 0, with 
7 > a, (3 > 6, 7 — a = (3 — 6, 5 — 7 = 7 — a + and that Z — A; + l = 7 — ain G'^p.^^{2) > 0. 
Therefore property (13.421) is still valid, and (13.451) is satisfied. An analogous computation 
to (^M-<^M yields 

sup |s^(6-,Ct-)l = sup(7 - a, 1,7 - a + 1) 
sup 1 5^(5^+26-, 5^,x+2Ct-) I = sup(l,7-a + 2,7-a + 1) 
CSiCtXt) > G^;;,^,,,(2)(sup |5.(^J.+26-,5:,.+2a-)| - sup|..(6-,Ct-)l) > 

thus (13.401) is wrong. 

(ii) Assume that y = x + 1, C,o{x — 1) = 7, ^0(2; — 1) = a, and that in ^^'^.^^(1), 
(3.48) 0<(3-6<-f-a<l-k<{-f-a) + {(3-6) 

so that (12381) is not satisfied while (12341) is. We have by (131421) . (ICTl) . ([336]) 

S{^t-Xt-) = s^_i(6,Ct) = Sx-i(6-,Ct-) =7-« 

Sx+i(6,Ct) = s^+i{^t-Xt-) = (3 - S 

Sx(6-,a-) = 0; s,{^t,Ct) = l-k = S{^t,Ct) 
so that (12.611) fails at time t because of (13.481) . I 
Proof of Proposition I2.64L 

Let (a, P, 7, 6) G X'^ and k, I non-negative. Assuming ( I2.65I) - (I2.66I) . we show that if k, I do 
not satisfy (1238|) . then GJJ^.^^ = 0. 

If / > is:, (12361) implies that S^_^ > S!^;/ and thus S^^^ - A S^^ ,^ > T^^^. Therefore 
G'a!(,;,,s = by (^M (if ^ > 0) or by (^M) (if = 0). ' 

If A; > L, (12351) implies that S^^ ,5 > S^;^^ and thus S!^ 5 - S^^^ A S^^ ,^ > L^^^. Therefore 

G^ai7,^ = 0by (imi) (if/>0)orby (ESni) (if/ = 0). ' ' ' ' ^ 

4. Coupling and attractiveness on examples 



In this section, we apply results of Section [2] to examples, first the classical ones, SEP- 
ZRP-MP, then a stick process (StP), and finally a two species exclusion model (S2EP). 
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4.1. SEP-ZRP-MP. According to notations in Sections [T] and [21 let p{x,y) = p{y — x) 
be an irreducible translation invariant probability transition on Z'^; set the transition rates 
T^jsiy — x) equal to unless k = 1, and write r^^^(?/ — x) = p{y — x)b{a,P). We specify 
below (cf. (12. ip ) the necessary conditions on the rates and the proper reduction of the state 
space (when needed) which ensure that fl2.15l) is the generator of a well defined Markov 
process (see [TS], [I], [S]) with a one parameter family of invariant and translation invariant 
product probability measures, denoted by {yUp}p. The parameter p is the average particles' 
density per site. We assume that the transition rates satisfy condition (A) of Section [H 

(a) MP: The misanthropes process ([H|)- There, X = N, S = Z"'. We assume that the 
transition rates b{a,f3),a,l3 e N are Lipschitz functions w.r.t. the first variable, with a 
Lipschitz constant uniform in the second variable; the state space is reduced to 

(4.1) ^0 = {r] e X^;^a{x)r]{x) < +00} 

where a is a function which satisfies J2yesP(^^y)^(y) — ^'^{^) x E S, for some 

constant M (see |lj, whose construction for ZRP extends for MP to these conditions, 
weaker than in [H (1.7)]). 

We quote [H (2.3), (2.4)], which imply the existence of product probability measures 
{fip}p, with p e [0, +00): 

(Ml) 6(0, .) = and b{a, /5) > for a > and /3 > (irreducibility on X^) 

( {yz e Z'^, p{z) = p{-z)) or 6(a, /3) - 6(/3, a) = b{a, 0) - b{(3, 0) 
(M2) I b{a,P) _ b{a,0)b{l,(3) 

(b{/3 + l,a-l) 6(/? + 1,0)6(1, a -1) 

(b) ZRP: The zero-range process ([1]) is a misanthropes process for which the rates depend 
on the first coordinate only, that is, b{a,f5) = g{a) (for more general existence conditions 
in the totally asymmetric case, we refer to [6]). Condition (Ml) reduces to 5^(0) = and 

(4.2) g{a) > for all a > 
Condition (M2) is verified in all cases and p G [0, lim^ (7(a)). 

(c) SEP: The simple exclusion process ([17], [18]) has state space {0,1}^''. The only 
nonzero jump rate is V\ Q^y — x) = p{y — x)6(l, 0) = piy — x). Irreducibility for p(x, y) is 
relaxed to 

(4.3) Va;, y G Z\ Wt > 0, pt{x, y) + pt{y, x) > 0, 

where pt{x,y) are the transition probabilities for the continuous time Markov chain as- 
sociated to p{x,y). The invariant product probability measures are Bernoulli, with p = 

pM^) = i)e[0,i]. 

In those three examples, (I2.22l) - fl2.23l) reduce to Condition (A), which thus insures 
attractiveness. For ZRP, (A) reads: g{a) is a non-decreasing function of a. For SEP, (A) 
is empty and does not impose anything. 
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The coupling rates' formulas fl2.27l) - fl2.28l) reduce to those of basic couphng, since the 
only non-zero possibility for k and Z is 1: 

^l;0 _ pi _ pi .pi 

In Section [6] we revisit the conditions previously stated to determine extremal invariant 
and translation invariant probability measures for SEP-ZRP-MP. 

4.2. A stick process (StP). This example is motivated by its connections with a gen- 
eralization of the discrete Hammersley-Aldous-Diaconis (HAD) process, a continuous-time 
Markov process taking values in a subset of A" = {0, 1}^ (see [IQ]). Let p{x, y) = p{y — x) 
be an irreducible nearest-neighbor translation invariant probability transition on Z, with 
p{l) +p{—l) = 1. For each empty site j, the closest particle to its left (resp. right) jumps 
to site j with rate p(l) (resp. p(— 1)). The formal generator of the generalized HAD process 
is 

(4.5) Lnfiv) = E E P(')[f(SUv) - fiv)] 

jez z=±i 

where, for all G A", x G Z and z = ±1, 

(4.6) Xz = Xz{rj,x) = zmax{y < zx : rj{zy) = 1} 

The existence of this dynamics is proved through a graphical construction, and to exclude 
the possibility of particles escaping immediately to infinity, the state space is restricted to 

{—1 m " 

rj E X : lim m~^^'^ rj{j) = lim ■m~^^'^ / ^viJ) = °o 
m-^+oo ' ^ m— >+oo ' ^ 

j=-m j=l 

Notice that particles cannot jump over each other, so they keep their relative order. We 
denote by {zi{t),i G Z, t > 0) their positions, with zo{0) the position of the first ?7-particle 
initially to the right of (or at) the origin. 

We study here a generalized stick process (StP) which is a discrete version with a gen- 
eralization to nearest neighbor jumps of the stick process defined in [20j (for a totally 
asymmetric case). It is defined on the state space 

-1 n 

fig = {?7 G N^; lim n~'^y^r]{x)= lim 71""^^^ r]{x) = 0} 

n— ►—oo ' ^ n— >+oo ' ^ 

x=n x=l 

SO that Qo C with A = N, 5 = Z; it describes translation invariant nearest-neighbor 
jumps. Its transition rates are 



(4.7) 

They satisfy (I2A8|) . 



ra,/3(^) = p{z)l{k<a} for z = ±1, 
'^a,t3i^) — otherwise 
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Similarly to the original models (see [20]), the generalized stick process is in bijection 
with the generalized discrete HAD process seen from a tagged particle: If we assume that 
at time 0, t/o(0) = 1 so that Zo{0) = 0, we have for all t > 0, z G Z, 

(4.8) r]t{t) = Zi+^{t) - z,{t) 

In other words, holes (resp. positions of particles) of the HAD are turned into particles 
(resp. sites) for the stick model. 

The geometric product probability measures {/ip, p G [0, +oo)} with parameter p(l+p)^^ 
are invariant for the generalized stick process (p is the average particles' density per site). 

Proposition 4.9. The stick process is attractive. The rates of the coupling generator C 
are, for all (a,/5, 7,5) G N"^, positive k, I, z = ±1, 

(4-10) G'ai7,5(^) = P{z)'^{k<a,l<J,a-J=k-l} 

(4.11) <^ai;7,5(^) = p{z)l{i<^-a} where 7 > a 

(4.12) G'a,°;7,5(^) = p{z)l{k<a-'r} whcrc a> 
This increasing coupling does not allow exchanges of discrepancies. 

Proof. Inequalities fl2.22p - fl2.23p are satisfied for the rates (14. 7p of StP, which is thus attrac- 
tive. The coupling rates formulas fl2.27p - fl2.29p write here fl4.10p - fl4.12p . Finally, recalling 
Theorem I2.52[ notice that the case A(a, /5, 7, 5, fc, /) = 0, (a,/?), (7,5) not ordered (line 2 
of fl2.55p ) is not possible; more generally (cf Definition 12. 56p exchanges of discrepancies are 
impossible since \k — l\ < \a — and thus a — j + I — k cannot change sign. I 

Remark 4.13. The rates fl4.10p - fl4.12p are very different from the basic coupling rates. 
For k, I positive, k = 1 in ^.^ ^{z) is possible if and only if a = 'y, and 

4.3. The two species exclusion model (S2EP). This one dimensional model with 
charge conservation was introduced and studied in [9j (see also [22], [23] for a full account 
on conservation laws in this model). Its state space is { — 1,0, 1}^, with the interpretation 
that a site on Z is either empty or occupied by at most one charge, positive or negative. A 
transition consists in a simultaneous change of the charges on two nearest-neighbor sites, 
such that the total charge is conserved: in a transition, the value at a site can change 
by -|-1, —1, +2 or —2, therefore the admissible values for k are 1 and 2. Since in [9] 
the model is interpreted as a Solid-on-Solid interface, conventions there are different from 
the present ones, and the rates r^^^(l) (resp. r^^^(— 1)) here are equal to their F"^ (resp. 
F^^^). There are possibly ten different transition rates, namely: Tl_-^^{z), Tl_^{z), T\ _i{z), 
TIq{z), T\q{z), z = ±1. We do not consider the case Fqq(z) = T\ _i{z) = 0, where there 
is a second conserved quantity (the number of particles involved in each transition). We 
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assume, in order to avoid degeneracies, that 

rU(i) + rU(-i)>o 



(4.14) 



r^,o(i) + r^,o(-i)>o 



Under the condition 

r^,o(-i)n,-i(i)-r^,o(i)rU(-i) ^ 
n,o(i) + r^,o(-i) 

(4.15) Yl ' (rtiW + rUW-r},oW-rtiW) 

z=±l 

the process has a one-parameter family of stationary product probabihty measures {/ip, p G 
[—1, 1]}, where fii = 6i and = 6-i are the Dirac measures concentrating respectively 
on ?7 = 1 and t] = —1. See Section [7] below for more details. 

Proposition 4.16. 1) The two species exclusion model is attractive if and only if the rates 
satisfy, for z = ±1, 

(4.17) Tl_,iz) V Tl.iz) < Tl_,iz) A r},o(^) < rl_,iz) V r},o(z) < r},_i(z) + Tl_,iz) 

2) The increasing coupling for the two-species model does not allow exchanges of discrep- 
ancies if and only if, for z = ±1, 

(4.18) Tl_,{z) <Tl,{z) 

Proof. 1. Inequalities (14.171) are just a compact form of fl2.22p - fl2.23p for this model. 

2. By Definition 12.561 we have exchanges of discrepancies for (x,y) G 5*^, (a,/9, 7,5) G 
with (a — 7)(/9 — 5) < 0, non-negative k,l with \k — l\ > |q; — 7I V |/3 — 5\. Since 
|a - 7I V 1/3 - 5| > 1, the only possibility is \k - l\ = 2 hence (fc, /) G {(2, 0), (0, 2)}. Also 
y — X = z has to belong to { — 1, 1}. For a < -y, P > 6 the only possible positive coupling 
rate is 

Gll,,-,{z) = rl_,{z)-rl,{z) 

and for a > 7, /? < 5 it is 

G'?;°i;o,o(^) = rti(^)-rio(^) 

I 

The rates of the coupling process are here rather different from those of basic coupling. 
Their explicit form is derived in Appendix (Section [8]). 

Remark 4.19. Hydrodynamic limits of non attractive two species exclusion models have 
been studied in [T2| (where the dynamics has two conserved quantities), and in [TT], Sec- 
tion 6], where the rates are Tl_^{l) = 2 + a,Tl_-^{-l) = a,rj_i(l) = rlo(l) = 1 + 
a, ri__i(-l) = ri_o(-l) = a, Tl,{l) =f3 + a, Tl^,{-1) = a, Tl^,{l) = a, Tl_,(-1) = (3 + a 
(for suitable positive parameters f3,cr), which do not satisfy fl4.17p . 

We determine in Section [6] {TCiS)^ and in Section [7] the hydrodynamic limit for StP 
and attractive S2EP under Assumption fl4.18p . 
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5. Invariant measures: Determination of {lr\S)^. 

The main result in this Section, Theorem 15.131 is the determination of the extremal 
invariant and translation invariant probability measures of the process {r]t)t>Q with gener- 
ator (12.151) . assuming it is attractive. This result, which does not depend on the existence 
of product invariant probability measures, was already known in some cases (for instance 
SEP). The proof's skeleton is classical (and goes back to \XIl), but the derivation of its 
main step ( (15.141) below) in a rather general setting is new; it sheds some light on the 
relevant properties of the model, and will enable (see Section [6]) to relax the assumptions 
on the rates on the examples ZRP-MP. This proof is linked with Section 12.31 

5.1. Preliminaries and statement of Theorem 15.151 We first give irreducibility condi- 
tions for the coupled process {C,t, Ct)t>o of Proposition 12. 25^ which are required assumptions 
for Theorem 15.131 Their goal is to ascertain that if a coupled configuration contains a pair 
of discrepancies of opposite signs, there is a positive probability (for the coupled evolu- 
tion) that it evolves into a locally ordered coupled configuration. Whenever basic coupling 
can be used (A; < 1), such event can be described as a motion of discrepancies towards 
neighboring sites where they eventually merge; as a consequence, irreducibility conditions 
are more intuitive and can be stated directly in terms of the marginal rates. Here more 
involved events may occur, such as (partial or total) exchanges of discrepancies (see Def- 
inition 12.561) and irreducibility conditions need to be expressed in terms of the coupling 
rates. 

Definition 5.1. A pair of sites {x,y) G S"^ is an edge for F, and we write {x,y) G E-p, if 



k a,f3eX 



> 



Remark 5.3. Because the left hand side of (15.21) is symmetric and the rates are translation 
invariant, we have 

(x, y)eEr^ (y, x)eEr^ (0, y - x) E Er 

Remark 5.4. Due to the definitions (12. 27p - (12.291) of the coupling rates, the set of edges 
Ey for r coincides with the set Eq of edges for G, defined as 

k,l a,l3eX'y,SeX 

Definition 5.5. For any pair of sites {x,y) G S"^, the distance d{x,y) is the number of 
edges in the shortest path of sites p between x and y, 

d{x,y) = inf |n : 3p G P"(x, y)} 

V{x,y) = {p = {xo,- ■ ■ ,Xn) e S""^^ : xo = x,Xn = y,"^! <i <n,{xi^i,Xi) E Er} 

Definition 5.6. For any edge {x,y) E Er, a finite path of coupled transitions on {x,y) 
starting from {a' , f3' ,6') E and ending at (a", 7", 5") G X'^ is a finite sequence 
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{iar,Pr,lr,5r)}o<r<u C for u E N \ {0}, suck that (q;o,/3o, 70,^0 ) = {a',P',7',S'), 

Pu, lu, Su) = (tt", P", l", S"), and for all < r < u - 1, 

Otr — Or+l = Pr+\ ~ Pt ='■ 7r ~ 7r+l = ~ ='■ ^r; OTld 

either {K > 0, > 0, G^'^XirA^V - ^) > ^) 
or {kr <0,lr< 0, G-^-Jl^^,{x - y) > 0) 

Remark 5.7. This notion is the natural 'building block' to construct an event in which 
the width of discrepancies decreases. It plays the same role as a single coupled jump in the 
basic coupling construction. 

We also define the set 

= {(a, (3, 7, 5) G : (a < 7, /3 > 6) or (a > 7, /3 < 5)} 

of values corresponding to couples of discrepancies of opposite signs (cf. Definition I2.46P 
located on two sites x, y of S, that is ^(x) = a,C,{y) = P, ({x) = 7, C{y) = 6. 

In order to cope with 'strongly' asymmetric processes (in the sense that r^.^(?/ — x) = 
for some {x,y) G Ej- and G X"^), we introduce the following notions: 

Definition 5.8. An oriented wedge is an ordered triple {xo,Xi,X2) G such that 
{xo, xi) G Er, {xi, X2) G Er and Xq ^ X2. 

Definition 5.9. Let Wr be a set of oriented wedges. For any pair of sites {x,y) G S"^, a 
PVr-path from x to y is a set of oriented wedges pwr = ("W^i) " " " ,Wn-i), n > 2, such that 
Wi = {xq,x\,X2) for all 1 < i < n — 1, x\ = x}^'^ , x\ = x]^'^ for all 1 < i < n — 2 and 
Xq = X, X2~^ = y. We denote by V^^{x, y) the set of all Wr-pathes with x, y as endpoints 
and n — 1 wedges. 

Definition 5.10. The set of sites S is Wr-connected if {S,Er) is a connected graph and 
for all {x, y) G S"^ such that d{x, y) > 2, there exists a W^-path from x to y or from y to x. 

Definition 5.11. If S is Wr- connected, for all {x,y) G S'^, we define the quantity dwr{x,y) 
as 

dw {x y) = l^^^'^"^ < 1 

^ ' I inf |n : 3p G y)} otherwise 

Remark 5.12. d^^ is a pseudo-distance since triangular inequality may not hold. If 
Wt is taken as the set of all oriented wedges, d{-,-) and dwr{-r) coincide. In general, 
we have only dwr{x,y) > d{x,y) since the existence of a Wr-path with n — 1 wedges, 
PWy = (^1) ■ ■ ■ ,Wn-i), Wi = {xl,x{,X2) for all 1 < i < n — 1, implies the existence of a 
path of sites p = (xg, x\,--- , x^~^, X2~^) with n -\- 1 elements and the same endpoints. See 
in Examples \6.6\ cases where dwri^^u) > d{x,y). 

The irreducibility condition (IC) for the coupled process is: 
(o) {S, Er) forms a connected lattice. 
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For all (x, y) G Ep, 

(a) for any couple of discrepancies {a, P, 7, S) G X|,, there is a finite path of coupled 
transitions on (x, y) starting from (a, /3, 7, 5) aioiig which f^^ decreases; 

(6) let Xii{y — x) = Xl{x — y) he the set of values e G X such that for any 
discrepancy with values {a, 7) located at x, there is a finite path of coupled transitions on 
{x, y) starting from {a, e, 7, e) which ends with a discrepancy on y; then either Xii{y — x) = 
X orXL{y-x) =X. 

When (b) is not satisfied, we replace (6) by {h') : 

There exists a set of oriented wedges Wr such that 

(6q) S is Wr-connected; 

for aU w = {xq, xi, X2) G Wr, Xr{xi - xq) [JXr^xi - X2) = X; 

(62) for all w = {xq,xi,X2) G Wr, all ei ^ Xr{xi — xo) and all 62 G Xr{x2 — xi), 
there is a finite path of coupled transitions on {xi,X2) from {€1,62,61,62) to (£3, £4, £3, £4) 
such that 63 G Xr{xi — Xq). 

Theorem 5.13. If the process {rit)t>o with generator (12.151) is attractive and satisfies As- 
sumption (IC), then 

1) if the state space of {rit)t>Q is compact, that is X = {a, - ■ ■ ,b} for {a,b) G Z^, we 
have (X n iS)g = {/ip, p G 71}, where 71 is a closed subset of [a, b] containing {a, b}, and fip 
is a translation invariant probability measure on Q with fip[r]{0)] = p; the measures Pp are 
stochastically ordered, that is, Pp < Pp> if p < p' ; 

2) if {r]t)t>o possesses a one parameter family {pp}p of product invariant and translation 
invariant probability measures, we have (Tfl S)^ = {pp}p- 

5.2. Proof of Theorem 15.131 We consider case (2). For (1), we refer to 0, Proposition 
3.1], which deals with a misanthrope type model The only difference to cope with here is 
fl5.14p below, the difficult part common to both cases. 

• The method. 

Let /i G (Xn and V be an extremal invariant and translation invariant probability 
measure for the coupled process {^t, Ct)t>o with marginals p and pp, for some value p of the 
parameter (see [H Lemma 4.4] or fiSi Proposition VIII.2.14]). We claim that 

(5.14) \f{x,y) G S', I7{(e,C) : ^x) > ({x) and ^y) < ({y)} = 

Then V{{^,C) : ^ < C or ^ > C} = 1 (see P Corollary 4.8], [HI Lemma VIII.3.2(b)]), 
so that p = ppg, where po = sup{p, p > pp} (see [H Proposition 5.5], [HI Proposition 
VIII.2.13 and Theorem VIII.3.9]), which concludes the proof. 

We are left to prove (15. lip , that is, there cannot be discrepancies of opposite signs under 
the coupled measure V, which we can also write 

(5.15) y{x,y)eS^, V(a,/3,7,5) G with a > 7,/3 < 5, 

m^, ■■ a^) = «, c(^) = 7, ay) = ay) = s} = o 
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For this, we show that in a coupled evolution where the initial configuration has at least 
on a pair of sites {x,y) G 5*^ a couple of discrepancies of opposite signs (a,/3,7,(5) G Xfy, 
then with a positive probability, the width of the positive discrepancy decreases to in a 
(succession of) coupled transition(s). The latter exists by Assumption (IC). 

We prove 05.151) by induction on the distance d{x,y) = n (resp. dw^{x,y) = n under 
Assumptions (IC,6')) between the discrepancies of opposite signs. 

To clarify the presentation, we turn some claims into lemmas, whose proofs are postponed 
to next subsection. 

• The case d{x,y) = 1. 

By Theorem \2.52\ the function f^y = ipx + i'y with = [^(-2) — defined 

in (12.481) satisfies (12.531) . that is £^x,yfxy{^iC) ^ 0- Since V is invariant and translation 
invariant for the coupled process, 

= / CipQdj7= / YCu,o'ipodiy + / Cp^yip^dly 

ves v£S 



Since each term is non-positive by (12.551) . J CQ^yf^^dV = for any site v, and by translation 
invariance we also get / Cx^yf^^ydV = 0. Using (13.391) . we obtain 

= v{ax)=a,ax) = l,m = PX{y) = S}x 
(5.16) Y (G^k.Ay - ^) + Gi':a;5,,i^ - ?/)) A(«, P, 7, S, k, I) 

fc-i^{0,(a-7)+{(5-/3)} 

We conclude with 

Lemma 5.17. Under Assumption (IC,a), equation (15.161) implies (15.150 for d{x,y) = 1. 
• The first steps to n > 2. 

We suppose (I5.15P satisfied for d{x,y) < n — 1 (resp. for d]y^{x,y) < n — 1 under 
Assumptions (IC, b')). Consider any path of sites p = (xq, ■ ■ ■ , Xn) G ^"(x, y), (resp. Wr- 

path pv^r = {{xo,Xi,X2),{xi,X2,X3), - ■ ■ , {Xn-2, Xn-1, Xn)) G Vw^{x,y)). FoT all 1 < 2 < 

n — 1, by the induction hypothesis, 

Hi^, C) : ^{Xo) > ({Xo),^{Xn) < C{Xn),^{Xi) > ({Xi)} 

(5.18) + z7{(e,C) : ^xo) > axo),aXn) < ((a^n), ^(a;.) < ax^)} = 
Therefore if for (ei, ■ ■ ■ , e^-i) ^ we define the event 

(5.19) ^p(ei,--- ,e„_i) = {(^,C) : ^{xq) = a > ({xq) = -f,^{xn) = P < C{xn) = 6; 

(VI <n-l, e(x,) = C(xi) =£.)} 
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we now have 

■ ^{xo) = a> C{xo) = 7 , ^{xn) =I3< CM = 6} 

= ^ V{£p{ei, ■ ■ ■ 

(£i,--,£„_i)eX"-i 

We show that the left hand side is zero by proving the existence of a path of sites p G 
?/) such that for all (ei, ■ ■ ■ ,en-i) E X""-*^, 

(5.20) I7(^p(£l,--- =0 

We fix arbitrary (ei,--- ,en-i) G and p = (xq, ■ " " ,Xn) G V"'{x,y) (resp. pwr = 

((Xo,Xi,X2), {Xi,X2,Xs), ■ ■ ■ , (x„_2,X„_i,X„)) G 

Lemma 5.21. Equation (15.201) is satisfied for (ei, ■ ■ ■ ,£:„-i) G X"''^ if either (a) 6i G 
Xr{xi - xo), or (b) e„_i G Xr(x„_i - x„). 

If moreover either X/j(xi — xq) = X or Xji{xn-i — Xn) = X Lemma (5. 21( a) or Lemma 
15.21( b) implies that fl5.20p holds for any (^i, ■ ■ ■ ,en-i) G and we are done. Other- 

wise, we have both Xji{xi—Xo) ^ X and XR(a;„_i — x„) 7^ X. Some more work is required, 
depending on which assumption we work with: 

• Assumption (IC, b) is satisfied. 

Thus Xji{xo — Xi) = Xji{xn — Xn-i) = X, and we construct a new path of sites p' = 
(xq, ■ ■ ■ , x'^) where 



(5.22) 



Xq xq X, x^ Xji y 

X[ = Xq + y - Xn-l 

<_i = y + XQ-xi 

^x'j^i — x'j = — Xj for all 1 < j < n — 1 

By translation invariance (cf. Remark 15.31) . p' G V"'{x,y). Since Xji{x[ — x) = X^iy — 
Xn~i) = X and Xr{x'^_i — u) = Xr{x — xi) = X, both cases of Lemma [5.211 apply for the 
path p', and T'{Sp'{ei, ■ ■ • , En-i)) = for all (ei, ■ ■ ■ , En-i) G X"~^. 
• Assumptions (lC,b') are satisfied. 

We consider the path of wedges pwr = {{xo,xi,X2), {xi,X2,xs), • ■ ■ , {xn-2,Xn-i,Xn)) G 
Vy{r^{x, y) and the induced path of sites p = {xq, xi, ■ ■ ■ , x„) G Vn{x, y). Using (IC, 6'^), we 
have 

(5.23) Xr{xi - Xo) U Xr{xi - X2) = X 
When n = 2, we conclude by Lemma [5. 2 1[ 

When n > 3, we assume that 

(5.24) £1 ^ Xr{xi - Xo), Sn-l i XR{Xn-l " x„) 

(otherwise Lemma [5.211 implies fl5.20p ). Then by fl5.23l) . 

(5.25) £1 ^ Xr{xi - Xo), Sn-l G XR{Xn-l - Xn-2) 
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Lemma 5.26. Under Assumption (10,62), Equation (15.201) is satisfied "^{si, ■■ ■ ,Sn-i) G 
such that ei ^ Xji{xi — Xq) and e^-i G X/j(a;„_i — Xn-'i)- 

By (15.251) . Lemma [5.261 exhausts the remaining cases for n > 3 and V{Ey,{e\^ ■ ■ ■ , Sn-i)) = 
for all (ei,-- - ,e„_i) G X"-i. 

Theorem 15.131 is proven. I 

5.3. Proofs of Lemmas. All three proofs share a common part that we state as another 
lemma. 

Lemma 5.27. Let {x,y) G -Er and a finite path of coupled transitions {{ar, f3r,'yr,Sr)}o<r<u 
on (x,y). Define a family of cylinder indicator functions {(pr}o<r<u on coupled configura- 
tions such that for all C) G fi^, 

rVO < r < n, M^, C) < X-tiOxlfiO 
Then J (pydV = for v = u — 1 or for v = u implies J ip^dV = 0. 

Proof. The key idea is: if there is a succession of coupled transitions from an event of 
I^-probability to another event, then the latter must also be of 77-probability 0. 

We make a downwards induction on < r < f. We start the induction with the 
hypothesis / (p^jdV = 0. Then we suppose that for < r < t> — 1 we have J (pr+idV = and 
write 



= y Ccjyr+idf 

x',y' a',f^' Y,S' 
k,l 

x',y' a',f3' Y,S' 

E G'^'pwAy' - y^, ^i'yC)rf^(^, C) 

k,l 

The integral over negative terms in the second line is zero as a function integrated on 
an event of I7-probability (by the induction hypothesis). The third line is a sum of 
nonnegative terms, so each of them is separately zero, in particular the term corresponding 
given in Definition 15.61 that is either 

kr > 0, Ir > 0,x' = X,y' = y, k = kr, I = Ir, «' = «r, P' = Pr, l' = 7r, ^' = ^r, 

G'^l-nrA^y-^) > 0' x~:/^(Oxi:f (c)0.+i(5fee,4yC) = M^,o 
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or 

kr < 0, Zr < 0, x' = y,y' = X, k = —kr, I = —Ir, a' = Pr, P' = C(r, 7' = Sr, S' = 'Jr, 

In both cases, this yields 

J (prdV = 

The induction is proven and Lemma [5.271 follows. I 
Proof of Lemma 15.171 

We have fixed an arbitrary element {a, f], 7, 6) of Xf^. By Assumption (IC, a), for any cou- 
pled configuration C) such that ^{x) = a, ^{y) = (3, C,{x) = 7, ({y) = 6, there is a finite 
path of coupled transitions {{ar, Pr,1r,Sr)}o<r<u on {x,y) starting from (ao, /^o, 70, (^o) = 
(a, P, 7, 6) such that 

(5.29) f^S^r^, S^-K) < ftyii^ 

where Ks = J2r=o ^r, Lg = ^^^g Ir for < s < u — 1. Without loss of generality, we assume 
u to be the first step of the path of coupled transitions in which /^^ decreases, and, by 
Theorem [2321 >^.'':,C) = /^(e, C) for all < r < m - 1, so that (a„ 7., Sr) e 

Xf), and equality (15.161) is valid for (a^, Pr, Ir, ^r)- For < r < m — 1, we use Lemma 15.271 
for 

By our assumption (15.291) on u and Definition 15. 6[ for r = m — 1, we have either 
fc, > 0, > and G'jjXnrA iv ' ^) > 0' 

A{ar, Pr, 7n Sr, kr, Ir) = 1^3^,1^, ^J^C) - /^.(^S"^' ^t'^ < 

or 

fc, < 0, < and G-^^;-^Jl^^ (x - y) > 0, 

A{Pr, OLr, Sr, OLr, —K, —U) = fXyi^xJy^i ^x^yC) ~ ft,y{3x,y~^^i ^x^y^^C) < 

In both cases, the sum in (I5.16P applied to (a.^, Pr,'yr,Sr) contains a negative term, thus is 
negative by Theorem I2.52[ Hence, (15.161) applied to {ar, Pr,1r,Sr) gives 

j ipu-idV = z7{(^, C) : ^(x) = au-i, C{x) = 1u-i, ^{y) = Pu-i, Civ) = Su-i} = 

For M > 1, by Lemma [5.271 with f = u — 1 we get 

J ^odV = z7{(e, C) : a^) = a, ((x) = 7, m = P, Civ) = S} = 

I 

Proof of Lemma 15.211 

We only derive Lemma 15.21( a). supposing that ei G Xii{xi — x) (The proof of Lemma 
15.21( b). when £„_i G Xn^Xn-i — y), is similar and omitted). Then, there is a finite path of 
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coupled transitions on {x,xi), {{ar, l3r,1r,Sr)}o<r<u such that (ao, /5o, 70, 5o) = {c(,£i,l,£i), 
Pr = for all 1 < r < M and > (by attractiveness {a,ei) > (7,£i) implies that 

For < r < M, we define the cylinder indicator functions 

= a^i) = V2 < z < n - 1); = p, ay) = s} 

These functions fulfill Conditions fl5.28p of Lemma 15.271 for the finite path of coupled 
transitions {{oir, Pr,'lr, Sr)}o<r<u on (x, Xi). In addition, we have 

< mx) ■■ a^i) < axi),ay) > ay)} 

Since d{xi,y) = n — 1 (resp. dwri^i^y) = ~ 1); 

ni^x)--a^i)<a^i),ay)>ay)} = o 

by the induction hypothesis (we assumed (15. 151) satisfied for d{x, y) < n~l under Assump- 
tion (IC, 6), resp. for dwri^, y) < n — 1 under Assumptions (IC, b')). Hence / ^p^'^dV = 0. 
Applying Lemma 15.271 with v = u implies then 



= y ipo'^'dv = v{£p{ei, ■ ■ ■ , €n-i)) 



Proof of Lemma 15.261 

We define 

jo = inf{j ■.l<j<n-l,eje Xnixj - xj^i)} 

By f l5.25p . jo exists and 1 < jo. Using (IC,b2), we make a downwards induction on 
1 < m < jo — 1 to construct a sequence of finite pathes of coupled transitions on {xm, Xm+i), 
\{tm ,fcm+i5^m 5 J jo<r<«„ oucn tnat 



and 



X 



-JO 



G Xr(x 



JO 



XjQ-l) 



e Xb{x 



starts the induction for m = jo — 1 (thanks to the definition of jo), and for each step 
1 < m < jo - 1 



(5.30) 



~(m,0) _ ~(m+l,-u„+i) 



and 

~(m,Um) 

Cm 



e 



m+l 
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For all 1 < m < jo ~ 1 ^"^^ all < r < Um, we define the cylinder indicator function 
^ri^X) = m,C)■■ax)=a,ax) = l,ia^^) = axi)=e,,yl<^<m-ly, 

iax^) = Cix.) = ef~''''^-^\ \/m + 2<t< jo); 

i^Xi) = ax^) = e,,Wjo + l<t<n-l); ^v) = P, Civ) = S} 

(with the expressions between parentheses possibly empty according to the values of m 
and jo) which verify fl5.28p and in addition 

(5.31) ^s^(e,c) = ^™+i(e,c) 

for all 1 < m < jo — 1. By (15.301) for m = 1, we can now use Lemma 15.21( a) to get 
J ^u-^dv = 0, then we apply Lemma (5.271 for v = ui which gives J ip^dV = 0. This is the 
first step (m = 1) of another induction to prove that J (p^dV = for all 1 < m < jo — 1. To 
go from m to m + 1, we use (15.311) then another application of Lemma [5.271 to {^™}o<r<Mm 
for V = Urn- For m = jo — 1, we get 

= y (p^^~^dV = V{£p{ei, ■ ■ ■ , En-i)) 

The Lemma is proven. ■ 

6. Irreducibility Condition and invariant measures on examples. 

In this Section, we determine what Assumption (IC) imposes on the rates of the various 
examples of Section HI so that Theorem 15.131 applies: For ZRP-MP, Assumption (IC) 
enables to extend the previously known results to probability transitions p{-, ■) such that 

(6.1) q{x, y) = satisfies (gSD 

For SEP, Assumption (IC) is proven equivalent to (14. 3p . Results on Stick process (StP) 
and on the two-species exclusion process (S2EP) are new. 

6.1. SEP. Theorem 15.131 for SEP is equivalent to [HI Theorem 3.9(a)]. Indeed 
Proposition 6.2. For SEP, Condition (14. 3 P forp{.,.) is equivalent to Assumption (IC). 

Proof. We assume Condition (14. 3p . and check whether Assumption (IC) is satisfied. The 
lattice being connected by (14.31) . (IC, o) holds; 

There exists two pairs of discrepancies of opposite signs 

X|, = {(0,1, 1,0), (1,0, 0,1)} 

By symmetry and Remark 12.331 we can restrict the analysis to the first element in X|, . By 
definition, for any edge {x,y) G Er, we have: 

< Go'li^oiy -x) + G'i;o°o,i(x -y) = p{y - x) + p{x - y) 
thus at least one of the two terms in the last sum is non zero. 
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• (IC, a). A path of coupled transition on (x, y) starting from (0, 1, 1, 0) along which f^y 
decreases is given by 

'{(0,1,1,0), (0,1,0,1)} ifp{y-x)>0 
{(0,1,1,0), (1,0,1,0)} ifp(x-2/)>0 

so that Assumption (IC, a) is always satisfied. 

• (IC, b), (IC, b'). The value of X^^y — x) can be read from (cf. ( I4.4p ) 

G 0,0,1 fliy -x) = Gilo (y -x) = p{y - x) 
Gi,li,o{x -y) = Gi'f i o(a; -y)= p{x - y) 



which gives 
(6.3) 



G XR{y — x) if and only if p{y — x) > 

1 G XR{y — x) if and only if p{x — y) > 

Thus Assumption (IC, b) requires 

(6.4) V(x, y) G Er, p{y — x) > and p{x — y) > 
If this fails, we take 

Wr = {{xo,Xi,X2) e : p{xi - Xo)pix2 - Xi) > 

(6.5) and either p{xi — X2) = or p{xo — xi)p{xi — X2) > 0} 

Thus (14.31) implies that for any pair {x,y) G S"^, there is a Wp-path from either x to y, or 
from y to X and (IC, fog) holds. 

Let {xo,Xi,X2) G Wr- By (16. 3p . G X^^xi — Xq), 1 G Xr(xi — X2) and (IC,6^) holds. 
Now suppose Xfi{xi — xq) 7^ X; then ^(xo — Xi) = 0, p(xi — X2) = and Xij(xi — xq) = 
-^i?(2;2 — Xi) = {0}. Since p(x2 — Xi) > 0, there is a finite path of coupled transitions from 
(1,0, 1,0) to (0, 1,0, 1) on (xi,X2) and (IC,62) holds. Assumption (IC) is thus satisfied. 

For the converse, we assume (IC). Assumption (IC, o) implies Condition (16. ip . By (16. 4p . 
(IC,6) implies then (14. 3p . If Assumptions (lC,b') hold for some Wr, let {xo,Xi,X2) G Wr 
be such that 1 ^ X/j( ). Assumption {lC,b[) implies 1 G Xptixi — X2), hence 

p{x2 — Xi) > thus G Xji{x2 — Xi) by (16. 3p : then Assumption (IC,62) implies that 
1 ^ Xr{x2 — Xi) (hence p{xi — X2) = by (16. 3p ). since there is no possible coupled 
transition starting from (1, 1, 1, 1). Similarly, if ^ Xr{xi — xq) (hence p{xi — xq) = by 
(El), by (IC, b[) we have G Xr{xi - X2) and p{xi - X2) > by (ED; then by (IC, b'^), 
^ Xji{x2 — xi) (hence p{x2 — xi) = by (16. 3p . thus 1 ^ Xr{xi — X2) and by (IC, b[) we 
have 1 G Xr(xi — xq) thus p{xi — Xo) > by (16. 3p ). since there is no possible transition 
starting from (0, 0, 0, 0). Therefore the only possible choice for Wr is (16. 5p . 

This implies (i^D. ■ 

Examples 6.6. Let S = 7? he endowed with its canonical basis (61,62). We denote by 
the origin. 
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1. Let p{ei) = Pi,p(-ei) = gi,p(e2) = P2,p{-e2) = Qi, with pi, qi,p2, q2 all distinct, 
Pi + Qi + P2 + (I2 = 1- There is an edge in between and its four neighoring sites hence 
Assumption (IC, 0) holds, and p{y — x)p{x — y)>0 for all {x, y) G Er so that Assumption 
(IC, b) is satisfied. 

2. Letp{2ei) =Pi,p(-ei) = gi,p(e2) =P2, with pi + gi +p2 = 1- 
Ev = {(0,2ei),(0,-ei),(0,e2)} 

W^r = {(0, -ei, ei), (0, 2ei,2ei + 62), (2ei,ei,0), (ei - e2,ei,0), (ei - 62, -62,0)} 

Notice first that the distances d and dw^ are distinct: for x = 0,y = 4ei — 62 we have 
d{x, y) = 3, dwri^, u) = 5, corresponding to the following path of sites and Wr-path 

P = (0,2ei,2ei - e2,4ei - 62) 
PWr = ((0, ei,2ei), (ei,2ei,2ei - 62), (2ei,2ei - e2,3ei - 62), 
(2ei - 62, 3ei - 62, 4ei - 62)) 

There is an edge in E-p between and its four neighoring sites hence Assumption (IC, o) 
holds. There exists a Wr-path between and each of its 4 second neighbors, Ci + 62, — ei + 
62,-61 —62,61 —62 so that (IC,6g) holds. By 06.31) . (IC,b[) is satisfied, and so is (10,62) 
(by inspection of all elements ofWr). 

6.2. ZRP. Applied to this case. Theorem 15.131 generalizes [H Theorem 1.9], since 
Proposition 6.7. For ZRP, Assumption (IC) is equivalent to (16.11) and (14. 2p . 

Proof. Because the transition rates have a product form and g{-) is non decreasing (by 
attractiveness), for all {x,y) G Ep necessarily p{x,y) +p{y,x) > 0, and for any 7 > the 
existence of a coupled transition on (x, y) from (0, 7, 7, 0) G Xj) reads 

{p{x,y) +p{y,x))g{-f) > 

that is Assumption (IC, a) implies (14.21) . the usual irreducibility condition on g{.). We now 
verify that (16. ip and (14.20 imply Assumption (IC). 

• (IC, a). Consider (x, y) G Er and a pair of discrepancies (a, /3, 7, 6) G Xfy on (x, y) such 
that a < 7, /9 > 5 (the other case being treated by reversing the roles of the two marginals). 
A path of coupled transitions on (x, y) starting from (a, P, 7, S) along which f^y decreases 
(by 1) is given by either 

(6.8) {(a,/3, 7, 5), ■ ■ ■ , (0,/5 + a,7 - a, 5 + a), (0,/5 + q;,7 - a - 1, 5 + a + 1)} 

when p{y — x) > (with rates p{y — x) times g{a), g{a — 1), ■ ■ ■ , g{l),g{'y — a)), or by 

(6.9) {(a,/?,7,(5),-- - ,(a + (5,/?-5,7 + 5,0),(a + 5+ l,/3-5 - 1,7 + 5,0)} 

when p{x — y) > (with rates p{x — y) times g{6),g{6 — 1), ■ ■ ■ , g{l),g{P — 6)). 

• (IC,6). Consider a < 7 and e > 0. Suppose p{y — x) > 0; then a path of coupled 
transitions on {x,y) starting from {a,e,'y,6) which creates a discrepancy on y is given by 

(6.10) {(a, e, 7, 5), ■ ■ ■ , (0, e + a, 7 — a, e + a), (0, e + a, 7 — a — 1, e + a + 1)} 
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(with rates p{y — x) times g{a), g{a — l), ■ ■ ■ , g{l), g{'j — a)), which leads to Xji{y — x) = X. 

I 

[H Lemma 4.6] constructs the path of coupled transitions fl6.8l) (resp. 06.91] ) to go from 
(a, P, 7, 6) to (0, /3 + a, 7 — a, 5 + a) when p{y — x) > (resp. to (a + 5, /5 — 5, 7 + 5, 0) 
when p{x — y) > 0). Lemma [5.211 corresponds to pj, Lemma 4.7]. 

6.3. MP. This is the third example with single particle jumps. Applied to this case, 
Theorem 15 . 1 31 generalizes P, Theorem 2.13]: we do not require product invariant probability 
measures, that is [8j (2.3), (2.4)], and we extend [El (1.7)]. 

Proposition 6.11. For MP, Assumption (IC) is equivalent to either (i) b{'y,a) > for 
all •y > 0, a > 0, and p{-,-) satisfies fl6.1l) .- or (ii) for all {x,y) G E^, p{x,y)p{y,x) > 
and &(■,■) may take some zero values along with 

(6.12) 6(7, a) >0 for ally >a>0 

In [8], (i) and (m) have to be simultaneously valid. 

Proof. As for ZRP, necessarily p{x,y) +p{y,x) > for all {x,y) G E-p. For any 7 > a the 
existence of a coupled transition on (x, y) from (a, 7, 7, a) G X|) reads (we cannot reduce 
ourselves here to a = since b{-, ■) depends on a, 7, and not only on a like g{.) for ZRP) 

{p{x,y) + p(?/,x))(6(a,7) A 6(7,a) + [6(7,a) - 6(a,7) A 6(7,a)]+ 

+ 7) - 7) A 6(7, a)]+) = (p(x, y) + p{y, x))b{y, a) > 

and thus (16.121) is valid. On the other hand, suppose that there is (x, y) G E^ such that 
p{x — y) = and that &(7o, 5o) = for some pair (70, 60) such that < 70 < Sq. By 
attractiveness, b{a,P) = for all (a,/3) such that a < 70 and jS > So and there is no path 
of coupled transitions on {x,y) starting from (a, /?, 70, (5o) G Xf). Hence Assumption (IC) 
excludes having simultaneously p{x — y) = for some {x, y) G E-p and nontrivial zeros in 
b{-, ■). We thus have the two (non mutually exclusive) cases (z) and (ii) stated above. 

Case (i) follows exactly the same lines as for ZRP. For case (ii), the construction of the 
correct pathes of coupled transitions on (x, y) G Ep depends on the values ao, (3o for which 
6(ao,/?o) = 0. To explain how to proceed, we show that a Misanthrope process with rates 
6(-, ■) such that 6(7, 5) > for all 7 > 5 and 6(7, 5) = for all 7 < 5 satisfies Assumption 
(IC). 

• (IC,a). Let (x, y) G Ep and (a,/9,7,5) G on (x, y) such that a < 7,/5 > 5. Suppose 
first that 7 > 5 and let r be the smallest integer such that 7 — r < 5 + r. Observe that 
j3 + r>5 + r>'-j — r>a — r and thus a — r+l<f3 + r — 1 while by definition of r, 
7 — r + l>(5 + r — 1. Now choose r' < r — 1 to be the smallest integer verifying both 
a — r' < 13 + r' and 7 — r' > 5 + r'; a path of coupled transitions on (x, y) starting from 
(a, /?, 7, 5) along which /^^ decreases (by 1) is given by 

{(a,/5,7,5),-- - ,(a-r',/? + r',7-r',(5 + r'),(a-r',/? + r',7-r'- l,5 + r' + 1)} 

(with rates j9(x, y) times 6(a, /?), 6(a — l,/5 + l),-- - , 6(a — r' + 1, /9 + r' — 1), b{y — r', 6 + r')). 
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Suppose now that 7 < 5. Take r as the smallest integer such that 7 + r > 5 — r. Hence 
r < (5 - 7 + l)/2 and since /5 > 5 > 7 > a, {(3 - r) - {a + r) > (5 - r) - (7 + r) + 2 > 1. 
Thus P — r > a + r and a path of coupled transitions on {x,y) starting from (a,/?, 7, 5) 
along which f^y decreases (by 1) is given by 

{(a,/3,7, 5), ■ ■ ■ , (a + r, /? - r, 7 + r, 5 - r), (a + r + 1,/? - r - 1,7 + r, (5 - r)} 
(with rates p{y, x) times b{6, 7), b{6 — 1, 7 + 1), ■ ■ ■ , b{6 — r + l,7 + r — 1), b{P — r,a + r)). 

• (IC, b). Consider a < 7 and e > 0. There are five possibilities and for each one we give 
a path of coupled transitions on (x, y) starting from (a, e, 7, e) which creates a discrepancy 
on y. 

1) e < a, ro = inf{r > : e + 2r > a} and 7 > e + 2ro : 

{ (a, e, 7, e), ■ ■ ■ , (a - tq, e + tq, 7 - tq, £ + tq), (a - tq, £ + tq, 7 - tq - 1, £ + tq + 1) } 

(with rates p{x, y) times 6(0;, e), b{a — 1, e+ 1), ■ ■ ■ , 6(a — ro + 1, e + ro — 1), 6(7 — ro, e + ro)). 

2) e < a and 7 = 5 + 2ro = a + 1 : 

{ (a, £, 7, e), ■ ■ ■ , (a - tq, e + tq, 7 - tq, e + tq), (a - tq + 1, £ + tq - 1, 7 - tq, £ + tq) } 

(with rates p(x, y) times 6(a, e), 6(a — 1, e + 1), ■ ■ ■ , 6(a — ro + 1, e + ro — 1), then p{y, x) 
times 6(e + ro, a — ro)). 

3) a < e < 7 : 

{(a, £,7,6:), (a, £,7- !,£ + 1)} 
(with rate p{x,y) times 6(7,6:)). 

4) 7 < Ti = inf{r > : e — 2r < 7} and a = e — 2ri = 7 — 1 : 

{{a,e,-f,e), - ■ ■ ,{a + ri,e - ri,-f + ri,e - ri),{a + ri,e - ri,j + ri - l,e - n + 1)} 

(with rates p(?/, x) times b{e, 7), 6(e — 1, 7 + 1), ■ ■ ■ , b{e — ri + 1, 7 + ri — 1), then p{x, y) 
times 6(7 + ri, e — ri)). 

5) 7 < e and a < £ — 2ri : 

{(a,£:,7,£), ■ ■ ■ , (a + ri,£: - ri,7 + ri,£: - ri), (a + ri + !,£ - ri - 1,7 + ri,£: - ri)} 

(with rates p(?/, x) times 6(e, 7), ^(e — 1, 7 + 1), ■ ■ ■ , b{e — ri + 1, 7 + ri — 1), b{e — ri, a + ri)). 
This leads to Xr{ii — x) = X. 

Assumptions (IC, a, b) are thus satisfied. ■ 

6.4. StP. For all {x,y) G Er (so that = 1), {a, 13) G with « 7^ 0,/3 ^ 0, by K7\ 

and the definition of £'r, 

5^ K,^(l) + rJ,J-l)) = p{l)a + p(-l)/? > 

k 

therefore necessarily p{l) + p{—l) > and Assumption (IC, o) holds. We now check that 
Assumptions (IC, a, b) are satisfied. 
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• (IC,a). A path of coupled transition on {x,y) starting from 7,5) G Xf^ with for 
instance a < •y, P > 6, along which decreases is given by (cf. f l4.10p - fl4.12l) ) 

■ { («, P, 7, S), {a, p,a,S + y-a)} if p(l) > (with rate ^"^."^(l) = P(l)) 
{ («, /?, 7, 5), (a + /? - 5, 6, 7, 6) } if p(-l) > (with rate ^^-'^^^.(-l) = p(-l)) 

• (IC, b). Consider a < 7 and e >0. Suppose p{y — x) = p(l) > 0; then a path of coupled 
transitions on (x, y) starting from (a, e, 7, 5) which creates a discrepancy on y is given by 

{(a,£:,7,e), («,£,«,£ + 7 - a)} (with rate G°'J.:^"(1) = p(l)) 

which leads to X/j(l) = X. When p{x — y) = p(— 1) > 0, a path of coupled transitions on 
(x, y) starting from (e, a, e, 7) which creates a discrepancy on x is given by 

{(e,a,£:,7), (e, a, a, e + 7 - a)} (with rate G^'J.":^" (-1) =p(-l)) 

which leads to Xl{1) = X. 

6.5. S2EP. We first give the general conditions on the rates under which Assumption (IC) 
is satisfied. They are obtained by direct inspection of all possible cases, we rely on Table 
[H To illustrate this, we present three examples which fulfill Assumption (IC) in different 
manners. Notice first that since the dynamics is nearest neighbor and non degenerated by 
f l4.14p . (IC,o) is always satisfied. 

• (IC,a) 

= {(-1,0, 0,-1), (-1,0, 1,-1), (-1,1, 0,0), (-1,1, 0,-1), (-1,1, 1,0), 

(-1, 1, 1, -1), (0, 0, 1, -1), (0, 1, 1, 0), (0, 1, 1, -1), (0, -1, -1, 0), (0, -1, -1, 1), 
(0,0, -1,1), (1,-1, -1,0), (1,-1, -1,1), (1,-1,0,0), (1,-1,0,1), (1,0, -1,1), 
(1,0,0,1)} 

The following expressions need to be positive for z = ±1. They correspond to coupled 
transitions on (0, z) from (a, /5, 7, 6) G X^ with a < 7, /5 > 5 for which decreases 

G°Ji',o;o,-i(^)+G'i\_i,o(-^) = rl-r{z)+rl-,{-z) 

1,-1 W+G^J^l; 1,-1 W + G};°,^-l,l(-^)+G'i;°i;-l,l(-^) 

= rl_,iz) + rl_,{z) + rl_,{-z) + rl_,{-z) 
GlL,oi^) + G\:%^,i-z) = r},^(^) + r},„(-^) 

G^^AWW +G'°:i,0;l,-l(^) +G'o:°i;-l,l(-^) = Tl-^{Z)+Tl_,{z)+Tl_,{-Z) 

GSi;w(^) + ,,_,(.) + G};°_,,,(-^) = Tl_,iz)+Tl_,iz)+Tl,i-z) 



We need also that either 

^0i0;l-l(^) + Gq0.i^_i{z) + Gqq._ii{-z) 

(6.13) = Tl_,{z) - [Tloi^) - Tl_,{z)]+ + rl.iz) A rl_,{z) + T^.i-z) 
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is positive for z = ±1 (this corresponds to a coupled transition on (0, z) from (0, 0, 1, —1) 
for which /q*"^ decreases) or 



(6.14) 



and (16.131) are simultaneously positive for 2; = 1 or for 2 = —1. There, f^^ decreases 
along paths of coupled transitions on (0, z) for which the first one induces an exchange of 
discrepancies (see Example I6.15l (m) below). 

. (IC,6), (IC,6') 

Up to an exchange of marginals, there are three possible values («,7) with a < 7 
for a discrepancy on a given site: (—1,0), (—1,1), (0,1). Therefore for each element 
e e {—1, 0, 1} and each site z e { — 1, 1} neighbor of the origin, a set of 3 inequalities has 
to be verified by the coupling rates whenever e belongs to Xr{z) 



-1 e x.(z) 



'G'°jtl;0,-lW>0 

i;l,-lW+G°il,-l;l,-lW>0 



or 

'Gl'',iz)>0 



y1;0 



yl;2 



. 1 '-^-1,0; 0,0 (^) + ^o',-l;0,o(~^) + ^o|-l; 0,0 (~^) > 



-yO;l 



1;0 



As previously, the 3 first lines correspond to coupled transitions which create a discrepancy 
on z, and the combination of the fourth and fifth lines corresponds to a finite path of 
coupled transitions on (0, z), {(0, —1, 1, —1), (—1, 0, 0, 0), (—1, 0, —1, 1)} which ends with a 
discrepancy on z. 



e Xb(z) 



^ G^_11q. q q{z) + Gq 'Li. q q{-z) + Gq \q q{-z) > 



■yO; 1 



1;0 



or 



G'o;'i;o,o(-^)>0 
I [Gl\,,-^iz) + Gl\.^,^_,{z) + Gl^'_,^^^^^^^ 

^-l,0;l,o(^) + ^o!-l;0,l(~^) > 
''^oio;l,o(^) + ^o!o;l,o(^) + ^0^0; 0,1 (~^) > 

or 

fG'J;o;i,oW>o 

I I G\:\ i,o(-^) + Gl\ i,o(-^) + Gl\, i,o(-^) + Gl\, i,o(-^) > 
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1 e Xn(z) 



or 



G};!i;i,o(-^)>o 

. [^o!o;l,o(^) + ^o',0;l,o(^) + ^o',0; 0,1 (~^) > 

Gt-i;i,ii-^)+Gl'_,,,,,{-z)>0 

Using the coupling rates' expressions given in Table [1] as well as attractiveness conditions 
fl4.17p leads to following conditions on the rates: 



-1 e Xr{z 



e Xn(z) 



1 G Xniz) 



or 

I [Tl^_,{z) - r?,_,(^) > and Tl,{z) + (T^.^i-z) - Tloi-^)) > 

f fTl^(z) + iTl_,i-z)-Tl,i-z))>0 
or 

ri,o(-^) > and iTl_,{z) + Tl_,{z) - Tl_,iz)) + Tl_,iz) > 

rUz) + Tl_,i-z)>o 
(rj^o(^)-ri,oW) + rj,o(-^)>o 

or 

[ [Tl.iz) > and (Tl^^i-z) + - ri,o(-^)) + r?,_i(-^) > 

f f n^ii-^) + rti(-^) - rl,o(-^)) + rti(-^) > o 

or 

n,oi-^) - rl-i(-^) > and iVl.iz) - Tl.iz)) + TU~^) > 
Ti_^(-z) + Tl_,i-z)>0 

[ri,{~z)>o 



Examples 6.15. We give three examples of coupling rates which satisfy fl4.17p (that is 
correspond to an attractive process), and Assumption (IC). 

(i) All the rates T'^^{z), z = ±1 are positive. Hence (IC,a) holds and Xfi{z) = X for 
z = ±1, so that (IC,6) is satisfied. 

(ii) All rates r^^(— 1) are equal to and all rates r^^(l) are positive (totally asym- 
metric case). Here (IC,a) holds only through fl6.13p for z = ±1 ifTlQ^l) < r^_^(l), 
and through the combination of fl6.13p . (I6.14p if < < rQQ(l). We always have 
-1 G XrIi),-1 i Xr{-\) and 1 G Xr{-\), \ ^ XrIi), therefore (IC,6) never holds. 
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Here Wr = {{x,x + l,x + 2) e Z^}. To have (IC, b[), we need in addition 

rri^o(i)>ri,o(i)^x«(i)D{o} 

< or 

[Tl^,il)>Tl,{l)^XR{-l)D{0} 

In both cases, (IC, 62) is satisfied. 
(iii) We choose: 

= r},_i(i) = Tl.ii), < Tl_,{i) = Tl_,ii) = rl,o(i) 
= r?,_i(-i), < rj,o(-i) = rj,-i(-i) = rl,o(-i) = rU(-i) 

i/ere many coupling rates are equal to zero, and so is (16.131) for z = —1; we have (IC,a) 
satisfied through the combination of (I6.13p . (16.141) for z = 1: the pair of discrepancies 
(0, 0, 1, —1) on (x, X + 1) can disappear only through paths of coupled transitions involving 
an exchange of discrepancies: 

{(0,0, 1,-1), (-1,1, 0,0), (-1,1, -1,1)}, or {(0,0, 1,-1), (0,0, -1,1), (-1,1, -1,1)} 

By inspecting all possibilities, (IC, b) is satisfied. 

For this model we can say more on invariant measures: 

Proposition 6.16. For the attractive two species exclusion model, under Assumption (IC), 
fip G le (for p G [—1, 1] under condition fl4.15p . or p eTZ otherwise). 

Proof. It is similar to the one of [H Proposition 1.1], but, as previously to derive (15.141) . 
we have moreover to use Assumption (IC). I 

7. HYDRODYNAMIC limits for the stick process AND THE TWO SPECIES 

EXCLUSION MODEL. 

In this section, we derive almost-sure hydrodynamics for these two nearest-neighbor 
one-dimensional models under Euler scaling following the constructive method of [4j (see 
also [2], [3], where hydrodynamics is proved in the usual sense), which does not require 
product invariant probability measures for the dynamics. We quote the result from [1] 
(after explaining a few preliminaries), then indicate the scheme of its derivation. The fact 
that more than one particle can jump at a given time induces only one modification to the 
original proof. 

7.1. The hydrodynamic limit result. A graphical construction of the models is needed; 
it is explained as follows. Since for both models the rates are bounded, the construction, 
restricted to compact state spaces in can be extended here. For each {x,z) G Z^, let 
{T^'^, n > 1} be the arrival times of mutually independent rate sup^gj^.^^^gj^^ ^{z) Poisson 
processes, let {U^'^,n > 1} be mutually independent (and independent of the Poisson 
processes) random variables, uniform on [0,1]. We denote by (X,JF, P) the probability 
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space corresponding to these families of variables. At time t = T^'^, the configuration rjt- 



becomes S'^^^^r]t- if 



and stays unchanged otherwise. Let also (Xq, J?-o, IPo) denote an 'initial' probability space 
large enough to construct random initial configurations rj^ = ?7o('^o) for ^ Xq. 

Let iV G N be the scaling parameter for the hydrodynamic limit. The empirical measure 
of a configuration t] viewed on scale is given by 

a^{7]){dx) = N~^y^^ri{y)5y/N{dx) 

It belongs to the set of positive, locally finite measures on M, equipped with the topology 
of vague convergence. 

Theorem 7.1. ^ Theorem 2.1] Let ■, N ^ N) be a sequence of Q-valued random 
variables on Xq. Assume there exists a measurable bounded profile Uo{.) on M such that 

(7.2) hm a^{ri^){dx) =Uo{.)dx, Po-a.s. 

Af— >oo 

that is, 



lim / tjj{x)a (?7q ){dx) = / ip{x)uQ{x)dx, Pq-o-s. 
Jr J 

for every continuous function ip on M. with compact support. Let {x,t) ^— > u{x,t) denote 

the unique entropy solution to the scalar conservation law 

(7.3) dtU + d,[G{u)]=0 

with initial condition Uq, where G is a Lipschitz- continuous flux function (defined in (17. 6p 
below) determined by the T's. Then, with JPq ® JP -probability one, the convergence 

(7.4) lim a^{T]Nt{Vo i^o)^^))id'X) = u{.,t)dx 

holds uniformly on all bounded time intervals. That is, for every continuous function ip on 
M with compact support, the convergence 

lim /v.(.)a."(,;y.)(<ix)= 

^^'^ Jr J 

holds uniformly on all bounded time intervals. 

For a; G Z, the microscopic current between sites x and a; + 1 is defined by 

k 

The macroscopic flux G is defined by 

(7.6) G{p) = ^i,[JoM] 
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for p G 7^ (a set introduced in Theorem I5.13p . and by linear interpolation on the comple- 
ment of TZ (which is at most a countable union of disjoint open intervals). 

The first step to prove Theorem 17.11 works here unchanged: it consists in first giving a 
variational formula for the solution of the hydrodynamic equation fl7.3p under Riemann 
initial condition 

(7.7) Mo(x) = Al{a;<0} + pl{x>o} 

(where A,p belong to TV), and then deriving hydrodynamics under this initial profile. 

The second step consists in extending the hydrodynamic limit to any initial profile uq 
through an approximation scheme. Its key tools are a finite propagation property both at 
microscopic and macroscopic levels (which work unchanged: for particles, we refer to [21 
Lemma 3.1] which does not rely on attractiveness) and the macroscopic stability property. 
The latter requires here a different proof, that we derived in Proposition 12.601 

In the following subsections, we compute the involved quantities, that is microscopic 
current, macroscopic flux, and check macroscopic stability for both examples. We saw 
that the stick process possesses a one-parameter family of product probabihty measures, 
and that for the two species exclusion it is the case under Condition (14.151] , otherwise the 
invariant measures are not explicit (cf. Theorem 15. 13p . 

7.2. Hydrodynamic limits for the stick process. For a; G Z, setting p(l) = p and 
p(— 1) = g, the microscopic current between x and x + 1 is 



Jx,x+l{v) = ^(pl{fc<>?(a:)} — 9l{fc<r;(x+l)}) 



= - r]{x){r]{x) + 1) - - r/(x + 1){t]{x + 1) + 1) 
so that the macroscopic flux 

(7.8) G(p) = ^p(l + p) 

is a convex function for p — q > 0. 

Proposition 12.601 holds here by Proposition 14.91 

7.3. Hydrodynamic limits for the two species exclusion model. Under conditions 
fl4.14p - fl4.15p . the model admits a one-parameter family of invariant product probability 
measures {pp}pg(_i^i) (cf. Section 14. 3p : the marginals of Pp can be written as (see [9], 
Equations (3.6)-(3.8)]) 

(7.9) ,M^)^1) - ^<'"+'' 



(7.10) ,.,(,,W=0) = — -J— — J. = 1 - 4,(p) 

1 + cy + cy ^ 



(7.11) pp(r/(x) 



1 + cy + cy~^ 2 
zy + cy 



cy + cy- 
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where y G (0, +00) is an auxiliary parameter related both to the mean charge, p 
J ri{x)dfip{ri) and to the mean squared charge, ip{p) = J ri{x)^dpp{ri) as 

ciy-y~^) , . c{y + y-'^) 



(7.12) p= '" 2 -1 ' "^^P^ 1^ u_ -1 

1 + cy + cy ^ 1 + cy + cy 

The constant c depends on the jump rates 



After elimination of y in (17.121) . (f{p) appears as a function of p G (— 1, 1) 



(7.14) ^(p) 



"^1^- 4c^2^' ^'^^ "^^^^ ^ V4c2 + p2(l -4c2) if 4c2 ^ 1 

if 4c^ = 1 



2 

Independently of (I4.15p . according to the values of the jump parameters, the scaling is 
either diffusive or Euler, and the model is studied either as a gradient model, or as an 
attractive one. In both cases, the microscopic current between x and x + 1 reads 

J.,.+M = x°:°+i(^)(ro,o(i)-rio(-i)) 

+Xx°+i('?)rl,o(i) - Xx'i+i(^)ri,o(-i) 
(7.15) +x°';+i(^)ri_i(i) - xi!i(^)ri_i(-i) 

+xi;+i(r?)(rU(i) + 2r?^_i(i))-x;iii(^)(rti(-i) + 2r2_,(-i)) 



7.3.1. Hydrodynamic limits in the diffusive case. This case, for which we refer to 
corresponds to the values of the jump rates under which the macroscopic flux computed 
under Euler scaling in (I7.19P satisfies G{p) = 0. We only derive the conditions under which 
the model is gradient: 



Lemma 7.16. The two species exclusion model is gradient if the rates satisfy 

r^.o(-i) = ri,o(i); r},o(-i) = rl,o(i); ri,-i(-i) = r^-iii) 
rU(-i) + 2r? ,(-1) =ri (i) + 2r? ,(i)= ri._,(i) + r}n (i) 



(7.17) 



Proof. The model is gradient if there exists a function q such that 

Because the right hand side of (17. 5p depends only on ri{x),ri{x + 1), Jx,x+i{v) has to be of 
the form 

9i{Xxl+i{v) - xii+iiv)) + 92{x%li{v) - Xx!xliiv)) + 93{Xx^liiv) " Xx^xliiv)) 
with gi + g2 = 93- This leads to (17.171) . and 

(7.18) q{v) = -[r},o(i) + ri_,(i)]^ + [ri_,(i) - tI,{i)]^ 
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We have used 



l{r,(x)=0} = 1 - l{»?(x)=l} - l{r^{x)=-l}] l{»?(x)=l} = Vi^) ^ ! l{r^ix)=-l} = Vi^) 




rj{x) — 1 
2 



to write in (17.151) the indicator functions as polynomials in //(x), ri{x + 1). 



7.3.2. Hydrodynamic limit under Euler scaling in the attractive case. Using expressions 
(I7.15P for the microscopic current and (I7.9I) - (I7.11I) for the product measures, we obtain 
the macroscopic flux as a function of the density p, 



It is a nontrivial function, a priori neither concave nor convex in the parameter p. Recalling 
(17. 6p and (17.141) . its second derivative reads 



C, = 2(rio(l) - rio(-l)) + 2c\Tl_,{l) - Vl_,{-1)) + Ac\Tl_,{l) - Vl_,{-1)) 

C2 = (ri,o(i)-ri,o(-i)) + (ri_i(i)-ri_,(-i)) 
Cs = (ri,o(i)-ri,o(-i))-(ri_i(i)-ri_,(-i)) 



According to the values of the jump rates, the flux function is either concave, or exhibits 
one or two inflexion points. This yields increasing or decreasing shocks, rarefaction fans 
as well as contact discontinuities for the unique entropy solution of the hydrodynamic 
equation (17. 3p with initial condition (17.70 . 

Proposition 12.601 holds under Condition (I4.18p . 

Example. We detail what happens for thermal hath dynamics. These can be compared to 
the simulations described in [H], though a discrete time Markov chain was used there. The 
model is used to describe the dynamics of a one dimensional Solid-on-Solid interface, that 
is a collection of heights on Z and a configuration r] describes only the height differences 
between neighboring sites. Jump rates are thus chosen according to the following rules: 

• Detailed balance holds with respect to a formal Solid-on-Solid Hamiltonian. 



(7.19) 



G{p) 



(rio(i)-rio(-i))(i-¥'(p))' 
+^(r},o(i)-r},o(-i))(i-^(p))(v^(p)+p) 

- ri-i(-i))(i - v{p)mp) - P) 
+7((rti(i) - rti(-i)) + 2(rti(i) - Tl_,{-i))){^{pf - p') 



(7.20) 




where 



6Xp{-f/'(S'^,j.+l77) — H{r])j — r^(^^^)^;.^^(^)_j, ( — 1) 
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for all ?7, all a; and all /c > 0, with 

There are two parameters: J > weights the height differences between neighboring sites, 
and E' > is an external field which drives the interface. 

• Rates are chosen invariant under the symmetry x — > —x, rj — >■ —rj. 

• Condition (14.15^ for the existence of invariant product probability measures holds. 

The above three conditions allow to write all jump rates in terms of two physically 
significant parameters a = exp(— i?/2) and b = exp(— J): 



rl,o(-i) = r^,-i(-i) = « 
i) = rti(i 
-i) = rU(-i, 

lrU(i) = rU(-i) = ^" 



0,0 
0,0 



V2 



V2 



a + a 



f7.211 



with 



so that 
(7.22) 



The macroscopic flux and its second derivative read 



+ a 



+ a 



-2 



k-2 



G{p) 
G"{p) 



,-2 



1 -4c2 
2c2 



which yields no inflexion point if c G [1/4, l/-\/2]. When c > l/-\/2 there are two inflexion 
points near p = ±1; when c < 1/4 the two inflexion points are near p = 0. The situation is 
thus rather similar to that of the discrete time model studied heuristically in [9]. We have 
derived rigorously all the cases they got using finite size scaling analysis. The situation 
in Figure [1] was beyond the reach of their analysis: for c < 1/4 and initial condition (17. 7p 
with p = —A and A larger than the abscissa of the inflexion point, the system develops two 
shocks moving with opposite velocities and separated by a rarefaction fan. 
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-0.4 -0.2 0.2 0.4 



Figure 1. Starting from an initial condition Uo{x) = A(l{a;<o} — l{a;>o}), 
the hydro dynamic hmit may develop two shocks going in opposite directions 
at constant speeds and a rarefaction fan in-between. Here c = .145 and 
A = .303. 
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8. Appendix 

In this part, we detail the derivation of couphng rates for the attractive two species exclu- 
sion model of Subsection I4.3[ In order to illustrate the coupling construction of Proposition 
\2.25\ we re-obtain the rates' values by direct inspection, without using equations (12.271) - 
fl2.30p but relying on the following previously stated construction rules: (o) marginals have 
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to be recovered; (i) in coupled jumps, the same departure and arrival sites are taken for 
both processes; (ii) partial order on departure and arrival sites has to be preserved; (iii) 
discrepancies should not increase. In some instances, these rules do not fully determine 
the coupling rates but give a family of solutions; a consistent choice can be done by the 
additional rule: (iv) the non-zero coupling rates are located on a "staircase-shaped path" 
in the (/c, /)-quadrant (recall Proposition I2.35p . Rules (o - iv) single out the solution given 
by ^M-^M- 

By rule (i), we fix the same initial site x and final site y with z = y—x = ±1, consider only 
the local configuration {^{x),^{y),({x)X{y)) as an element of X^, with X = {—1,0,1}, 
and describe a coupled jump as a transition on this reduced state space. Furthermore, since 
in the initial process, jump rates from x to y depend only on values at these sites, one can 
look for a coupling process which has also this property so that the coupling construction 
can be given as if it were on X^. In Figure [2|, we represent this reduced space as an x 
square array with the values of {^{x),^{y)) (resp. (C(a;), Civ))) as the first, horizontal (resp. 
second, vertical), coordinate from left to right on the horizontal axis, and from bottom to 
top on the vertical one. Before asking for the actual values of the coupling rates, note that 
three different factors limit the number of non zero coupling rates: 

(a) there is charge conservation on each marginal. Since both ^(x) +^{y) and C(^) + C(?/) 
take value in {—2, ■ ■ ■ , +2}, charge conservation splits the x X^ square array into 25 
sub-arrays, and a transition takes place only between two elements of the same sub-array; 

(b) only positive charges jump from x to y in both marginals so that in a jump, by rule 
(i), the values at x cannot increase; 

(c) the coupling should be increasing, reflecting that the initial process is attractive by 
rule (ii) and that by rule (iii), the discrepancies cannot grow. Both properties are related 
to the function (cf. fICTD ) 

Once fixed the values of the local charges ^(x) + ^{y) and ({x) + ({y), fx,y{^X) takes its 
extremal values on ordered pairs, namely. 



if (e(x),e(y))<(C(a:),C(y)) 

e(x) + ^y) - (C(x) + ay)) if (e(x), e(y)) > (C(a:), ay)) 



Rule (ii) requires that f^y{^, C) cannot escape its minimum by a jump between x and y, 
and rule (iii) extends this requirement to non ordered states by asking that f^^yiC,, () is not 
increasing. 

These three restrictions can be depicted on the graphical construction of Figure [21 where 
we represent a possible coupled transition as an arrow from a 'cell' (initial values) to another 
one (final values). By (a), an arrow can be drawn only between two cells of a same subarray. 
By (b), since configurations are ordered by increasing charge and increasing value, arrows 
have to point left and/or downwards; In order to use (c), we have written in each cell the 
value of f^^y{^, ()■ arrows cannot go from a given value to a larger one. 
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All this already fixes trivial rates for our coupling process, when only one transition from 
a cell is allowed, and we have indicated all of them on Figure [H these include non coupled 
^-jumps (horizontal arrows), non coupled ^-jumps (vertical arrows) or fully coupled jumps 
between identical configurations (arrows on the main diagonal of each sub-array). We are 
thus left with only four types of coupled jumps, which we treat in more detail now, using 
the four examples represented in Figure [H 
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Figure 2. The coupling rates between two sites for the two species model. 
Possible coupled jumps, represented by arrows, are subject to three restric- 
tions: arrows connect two cells in the same sub-array, point leftwards and/or 
downwards, cannot point to a larger value than the initial one. All trivial 
arrows have been represented, leaving four non trivial cases. 



In the sub-array Ai = {(0, 1), (1, 0)} x {(—1, 0), (0, —1)}, all four pairs are ordered and 
the function fx,y{C, C) is constant. Starting from the state ((1, 0), (0, —1)), three jumps are 
possible and by rule (o), the coupling rates are related through 

^l!0;0,-l(^) = ~ G'i'o;0,-l(^) 
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Figure 3. Examples of non trivial cases 
Any choice for (j'i'o-o which preserves non-negativity is valid, so that 

Here the only possible non- zero value for k and Z is 1. Rule (iv) imposes that either 
o o-i(^) is zero, hence 

and the non-zero coupling rates are either on the (A;, /)-path Pi,o;0,-i = {(1; 0); (1) 1)} 
or Pi,o;o,-i = {(0, 1), (1, 1)}, depending on the value of the marginals. Notice that the 
construction is here identical to the 'basic coupling' one. 

In rectangle A2 = {(0, 1), (1,0)} x {(—1, 1), (0,0), (1, —1)}, we consider coupling rates 
from the initial state ((1,0), (0,0)). Its coordinates being ordered (that is (1,0) > (0,0)), 
order preservation (rule (ii)) forbids a jump to ((1,0), (—1, 1)), thus 

^i!o; 0,0(2^) = 

The only choice to get the correct (^-marginal (that is to satisfy rule (o)) is then 

^1^0; 0,0 (^) = ^0,oi^) 

and the last rate follows 

^l',0;0,o(^) = 0(2;) - Tl oiz) 

It has to be non- negative, which was stated by (14.171) . Notice that rule (iv) is automatically 
satisfied with Vifi-o^o = {(1, 0), (1, 1)}. 

The situation in rectangle A3 = {(-1, 1), (0, 0), (1, -1)} x {(-1, 0), (0, -1)} is simi- 
lar to the preceding one. We consider there the other type of non trivial initial state, 
((1, —1), (0, —1)) whose coordinates are again ordered ((1, —1) > (0, —1)), and by rule (ii) 

Gl\.^,^_,{z) = Gllliz) = 
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Here the situation cannot be solved through basic couphng, since there is a marginal ^-jump 
involving two positive charges. There is only one solution for the coupling rates, 

Gj;!i;o,-i(^) = Tl.^iz) + Tl_,iz)-Tl_,iz) 

They are all non-negative by (14.171) and identical to the values given by fl2.27l) - fl2.30l) . Rule 
(iv) is again satisfied and "Pi^-ijo.-i = {(1, 0), (1, 1), (2, 1)}. 
We find a last type of non trivial case in the larger square 

^4 = {(-1, 1), (0, 0), (1, -1)} X {(-1, 1), (0, 0), (1, -1)} 

and we consider ((1, — 1), (0, 0)) as initial state. Its coordinates are not ordered so that, 
rather than rule (ii), selection rule (iii) may apply. This does not fix completely the coupling 
rates and we have 

Gl;\o,oW = Tl-,iz) + Tl_,iz)-Tl,iz)-Gl\.^,^,iz) 

where non-negativity is guaranteed by (14.171) . and any choice for Gi'\.qq{z) in the range 

[Tl_,{z) - Tl,{zr < G?;°i^o,o(^) < rl_i(^) A {Tl_,{z) + Tl_,{z) - Tl,{z)) 
is valid. Rule (iv) imposes to set it to its minimal value, 

and gives a path Pi,-i;o,o = {(1, 0), (1, 1), (2, 1)} or Pi,-i;o,o = {(1, 0), (2, 0), (2, 1)} depend- 
ing on the value of the marginals. The solution given by (I2.27p - (l2.30p is again recovered. 
Note that when rl_'^^{z) < TqqIz), this choice is the only one which avoids an 'exchange 
of discrepancies' (cf. Proposition 14.161) . 

We conclude that (I2.27I) - (I2.30I) lead, as in the spirit of 'basic coupling', to changes as 
similar as possible in both marginals, which has in addition the property of minimizing 
exchanges of discrepancies. 

For the sake of completeness, all the coupling rates for this model are reported in Table 
[T| they are also useful to derive conditions on the rates in Subsection 16. 5[ 
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Table 1. non zero coupling rates for the two species exclusion mode! 



